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ABSTRACT 

Aims. The aim of this study is to enhance our knowledge of the characteristics and distribution of the circumstellar dust associated 
with the individual components of the young hierarchical triple system T Tau. 

Methods. To reach this goal, observations in the N-band (8-13 fim) with the two-telescope interferometric instrument MIDI at the 
VLTI were performed. For the northern component of the T Tau system, projected baseline lengths of 43 m, 62 m, and 85 m were 
used. For the southern binary projected baseline lengths of equivalent resolution could be utilised. Our study is based on both the 
interferometric and the spectrophotometric measurements and is supplemented by new visual and infrared photometry. Also, the 
phases were investigated to determine the dominating mid-infrared source in the close southern binary. The data were fit with the 
help of a sophisticated physical disc model. This model utilises the radiative transfer code MC3D that is based on the Monte-Carlo 
method. 

Results. Extended mid-infrared emission is found around all three components of the system. Simultaneous fits to the photometric 
and interferometric data confirm the picture of an almost face-on circumstellar disc around T Tau N. Towards this star, the silicate 
band is seen in emission. This emission feature is used to model the dust content of the circumstellar disc. Clear signs of dust 
processing are found. Towards T Tau S, the silicate band is seen in absorption. This absorption is strongly pronounced towards the 
infrared companion T Tau Sa as can be seen from the first individual N-band spectra for the two southern components. Our fits 
support the previous suggestion that an almost edge-on disc is present around T Tau Sa. This disc is thus misaligned with respect to 
the circumstellar disc around T Tau N. The interferometric data indicate that the disc around T Tau Sa is oriented in the north-south 
direction, which favours this source as launching site for the east-western jet. We further determine from the interferometric data the 
relative positions of the components of the southern binary in the N-band. We find good agreement with recent position measurements 
in the near-infrared. 

Key words. Stars: individual: T Tauri - Stars: circumstellar matter - Stars: pre-main sequence - Techniques: interferometric - 
Infrared: stars 



1. Introduction 

Our knowledge of the T Tau triple system has grown as obser- 
vational techniques were improved. T Tau was chosen by Joy 
(1945) as the prototype for a new class of variable objects. The 
members of this class are characterised as low-mass (< 3 M ) 
pre-main sequence stars still surrounded by accretion discs and 
thus exhibiting large infrared excesses (Appenzeller & Mundt 
1989). A number of nebulous patches have been detected in the 
vicinity of T Tau. They show similar spectral features to Herbig- 
Haro objects. Surrounding the source on scales of arcseconds, 
they clearly witness stellar outflows. Today an east-west and 
a northwest-southeast oriented outflow, as well as many inter- 
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locking loops and filaments of H2, are known, e.g., Herbst et al. 
(2007). 

A close companion was detected by Dyck et al. (1982). It 
was later named T Tau S because it is located 0.7" south of 
T Tau N. T Tau S is the source of one of the jets and a pro- 
totypical infrared companion (IRC) deeply embedded and not 
detectable in the visual (Stapelfeldt et al. 1998). Various expla- 
nations for its higher extinction have been discussed. On the one 
hand, a different evolutionary status of the southern component 
may be the reason. Already Dyck et al. (1982) saw in T Tau S 
a star that is less evolved than T Tau N. On the other hand, the 
extinction may be caused by a special geometry of the system 
like viewing T Tau S through an edge-on disc or viewing a back- 
ground southern component through a foreground screen. Then 
T Tau S can be a normal T Tauri star coeval with its northern 
neighbour. Koresko et al. (1997) concluded T Tau S is experi- 
encing episodes of embedding, perhaps because of an enhanced 
accretion rate at certain orbital phases. 
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T Tau S was resolved by Koresko (2000) into two compo- 
nents named T Tau Sa and T Tau Sb, with a separation of about 
0.05" in 1997. While T Tau Sb appears to be a "normal" ac- 
tive low-mass pre-main-sequence star simply residing behind 
an absorbing screen with Av ~ 15 mag (Duchene et al. 2005), 
T Tau Sa remains an enigmatic source. According to the picture 
drawn by Duchene et al. (2005), T Tau Sa is the most massive 
star of the T Tau system with a stellar mass of 2.5 - 3 M . A 
large fraction of the visual and near-infrared light emitted by the 
photosphere of this intermediate-mass young stellar object are 
scattered away thanks to an edge-on circumstellar disc or lost 
through its "polar" regions. Furthermore, a circumbinary enve- 
lope or thick disc may obscure both southern components lead- 
ing to an additional visual extinction of Ay ~ 15 mag as it has 
been found towards T Tau Sb. 

Associated with T Tau Sb is a radio source observed by 
Loinard et al. (2007) with the VLBA to determine its distance 
precisely. The derived value of 146.7 + 0.6 pc has been adopted 
in this paper for the whole T Tau system. 

The northern component has shown a significant brightness 
variation neither in the K- nor in the L'-band (Beck et al. 2004). 
Even its visual flux seems to be constant (Ghez et al. 1991). The 
southern component shows variability with a total amplitude of 
~ 2 mag in the K-, and ~ 3 mag in the L'-band. Beck et al. (2004) 
explained the variations by changes in the amount of obscuring 
material. In its bright state, T Tau S dominates the flux of the sys- 
tem longwards of 3 //m. According to Ghez et al. (1991), the flux 
ratio between the southern and the northern component in the 
N-band ranges from ~ 0.4 in the minimum to ~ 2.6 in the max- 
imum. Both components forming T Tau S have been identified 
as variable. The flux ratio of T Tau Sb with respect to T Tau Sa 
varies between ~ 0.1 (Koresko 2000) and ~ 3.6 (Beck et al. 
2004) in the K-band. 

There has been general agreement that T Tau N is surrounded 
by a circumstellar disc of moderate mass (Beckwith et al. 1990) 
seen not far from face-on, although its proposed size (Akeson 
et al. 1998) is still under discussion. For T Tau S it was not at all 
obvious that its circumstellar material would be concentrated in 
a circumstellar disc. Its outburst in 1989 led Ghez et al. (1991) to 
model its SED as dominated by a FUOR-type accretion disc. The 
picture of T Tau S then shifted to that of a rather normal young 
star (Hogerheijde et al. 1997) extincted by a disc seen nearly 
edge-on (Solf & Bohm 1999), but with very low mass according 
to the limits derived by millimetre interferometry (Hogerheijde 
et al. 1997; Akeson et al. 1998). This missing mass was natu- 
rally explained by truncation due to the companion T Tau Sb 
(Koresko 2000). Ultraviolet studies from HST showed an ab- 
sorbing screen in front of T Tau S (Walter et al. 2003), possi- 
bly associated with a circumbinary disc around T Tau S, which 
would be responsible for most of the extinction towards T Tau 
Sb. By mid-infrared adaptive optics imaging in the silicate fea- 
ture Skemer et al. (2008) add to the evidence for different disc 
orientations in the T Tau system and question the importance 
of a circumbinary disc. The actual orientations of the circum- 
stellar discs in the close binary T Tau S have remained largely 
unknown. 

Summarising the brief description above, the T Tau system is 
one of the best-suited laboratories for investigating the dynamics 
and the formation of young multiple systems. A first step in this 
direction has been made by monitoring the orbital motions of the 
southern components around each other and around the northern 
component. These studies led to reliable masses of the stars and 
concluded that T Tau Sb is indeed bound to the system (Duchene 
et al. 2006; Kohler et al. 2008; Kohler 2008). A next step is the 



investigation of the individual discs with respect to their dust 
content and their geometry. The information about the geometry 
will help in addressing the issue of coplanarity of the discs with 
respect to each other, as well as to the orbits. 

In this paper we present observations (Section 2) of the 
T Tau system with the "MID-infared Interferometric instru- 
ment" MIDI, which resides in the interferometric laboratory of 
the "Very Large Telescope Interferometer" (VLTI) on Paranal, 
Chile. After a brief description of the data reduction (Section 3), 
we present first the directly derived results (Section 4). In 
Section 5, we model the northern component with a radiative 
transfer model and investigate the dust properties of its circum- 
stellar disc by fitting the correlated fluxes. In Section 6, the vis- 
ibilities of T Tau S are analysed to derive relative positions and 
separate spectra for the two components. With the help of radia- 
tive transfer models, the properties and the relative orientations 
of the discs are determined. Finally, a picture of the system is 
drawn in Section 7. In Section 8 a summary can be found. 



2. Observations 

We obtained three separate interferometric measurements of 
both T Tau N and T Tau S within the "Guaranteed Time 
Observations" (GTO) in October and November 2004. The base- 
lines used were UT2-UT3 (Oct. 30/31), UT2-UT4 (Nov. 1/2), 
and UT3-UT4 (Nov. 3/4). A journal of observations appears in 
Table 1. 

2.1. Observing Sequence 

The standard observing sequence for an interferometric mea- 
surement with MIDI on the VLTI produces an image of the ob- 
ject at 8.7/zm, a spectrum from 8/zm to 13 /urn, and spectrally 
resolved visibilities and correlated fluxes over the same wave- 
length range. For descriptions of the instrument and its opera- 
tion see Leinert et al. (2003a,b), Morel et al. (2004), and Ratzka 
(2005). 

The chopped image is taken after the coarse acquisition by 
the telescopes. These images allow the observer to adjust the 
position of the objects to a predetermined pixel in order to max- 
imise the overlap of both images for the following interferomet- 
ric measurements. 

Then, the beam combiner, which produces two interferomet- 
ric outputs of opposite sign, is put into the optical train. Together 
with a dispersion device - in our case, the low resolution prism 
(A/AA « 30) - dispersed fringes, i.e., spectrally resolved inter- 
ferograms, can be obtained. After locating the zero optical path 
difference (OPD) by scanning a few millimetres around the ex- 
pected point of path length equalisation, an interferometric mea- 
surement with self-fringe-tracking is started. In this mode, the 
piezo-mounted mirrors within MIDI are used to scan a range in 
OPD of 40-80 fim in steps of typically 2 /im. At each step, with 
the corresponding fixed OPD, an exposure is taken that gives 
the instantaneous value of the interferometric signal from 8 pm 
to 13 yum. After each scan, the position of the fringe packet in 
the scan is measured and the VLTI delay lines are adjusted in 
order to recentre the fringe packet for the next scan. Even for 
the low resolution of the prism, the coherence length is about 
+ 300 /im. This means that imperfect centring of the interfero- 
metric scan within a few A has no adverse effect on the visibility 
determination, i.e., no noticeable spectral smearing occurs. The 
output after a scan is the spectrally resolved temporal fringe pat- 
tern, which gives the fringe amplitude or correlated flux. The 
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Date of 


Universal 


Object 


IRAS 


Proj. Baseline* 


Airmass** 


Acquisition 


Interferometric 


Photometric 


Observation 


Time 




[Jy] 


[m] 


[deg] 




Frames 


Frames 


Frames 


.31-1U-ZUU4 


03:53 


- 04:06 


riD ZU044 


i a n 
14. / 


OO 1 

zo. 1 


C/1 O 

54.Z 


1 oo 
l.oz 


1 HAHv/l mo 

1UUUX4 ms 




1 Cf|A\/ 1 O mo 

zxi juux iz ms 


.31-1U-ZUU4 


06:45 


- 06:59 


riD ZjoU4 


c i 
J.L 


/ion 
4Z.U 


5U.4 


i ^n 
l.DU 


1UUUX4 ms 


i nnnnv 1 o mo 
1UUUUX1Z ms 


Tvonn/iv 1 o mo 

zxzuuux 1 z ms 


q i in onn/i 

j 1-1U-ZUU4 


n£. ^n 
Uo. Jy 


m.oo 
- U/.ZZ 


T Ton M 

i iau in 




/ton 


a n a 
49.4 


1 /1< 
1.4 ' 


1UUUX4 ms 


1 (\C\(\(\^/ 1 mo 

iuuuux iz ms 


Ovonnriv 1 o mo 
zxzuuux i z ms 


1 1 in onn/i 
jl-lU-ZUU4 


07:22 


- 07:41 


T TV,, , C 

1 lau o 




43.9 


48.9 


i a n 
1.49 


1UUUX4 ms 


i nnnnv 1 o mo 
IUUUUX iz ms 


Tvonn^v 1 o mo 

zxzuuux 1 z ms 


31-10-2004 


07:41 


- 08:10 


HD 37160 


6.5 


43.5 


46.2 


1.22 


1000x4 ms 


10000x12 ms 


2x2000x12 ms 


31-10-2004 


08:40 


- 09:08 


HD 37160 


6.5 


45.8 


46.7 


1.31 


1000x4 ms 


10000x12 ms 


2x2000x12 ms 


no 1 i onn/i 

UZ-l 1-ZUU4 


05:10 


-05:36 


T Ton M 

i lau IN 




o c n 


q-7 n 
8/. 9 


1 /10 
1 .4Z 


1 nnnv/i mo 
1UUUX4 ms 


i ^r\r\r\\y 1 o mo 
louuux iz ms 


o \^ o ^ nn vii mo 
zxzjuux iz ms 


no 1 i onn/i 
UZ-l 1-ZUU4 


05:36 


-06:10 


1 lau o 




s /.o 


or /: 

85.0 


i /in 
1.4U 


1UUUX4 ms 


i Annn^ 1 o mo 
loUUUX IZ ms 




no 1 i onn/i 

UZ- 1 1-ZUU4 


06:10 


- 06:37 


| I p, Q 1 A 1 
1 1 \J 314Z1 


Q ^ 


0*7 1 


C/1 ^ 


i "in 
Y.LI 


1UUUX4 ms 


i /^nnriv 1 o mo 
louuux iz ms 


OvT^nriv 1 o mo 
zxzjuux iz ms 


no 1 i onn/i 

UZ- 1 1-ZUU4 


07:08 


- 07:41 


1 1 \J J I I OU 


O.J 


on 1 
89.1 


oo o 
8Z.S 


1 1 

1.Z1 


iuuuux4 ms 


i Annriv 1 o mo 
louuux iz ms 


o \^ o ^ nn vio mo 
zxzjuux iz ms 


no 1 i onn/i 
UZ-l 1-ZUU4 


08:09 


- 08:37 


tiU 3U / la 


1 1 .5 


88.6 


80.2 


i nQ 


1 flAHv/l mo 

1UUUX4 ms 


1 AAHAv 1 O mo 

loUUUX IZ ms 


Ovoc nn wio ,-v-, r, 

zxz juux iz ms 


no 1 i onn/i 
UZ-l 1-ZUU4 


08:37 


- 09:03 


i in a i no c 


n 1 
/.I 


87.2 


79.6 


i n^ 


1 flAHv/l mo 

1UUUX4 ms 


1 /;nnnv 1 o mo 
loUUUX IZ ms 


OwTc nn v/H mo 

zxz juux iz ms 


04-11-2004 


00:01 


- 00:26 


HD 178345 


8.6 


57.0 


145.6 


1.42 


1000x4 ms 


20000x12 ms 


2x5000x12 ms 


04-11-2004 


02:10 


- 02:47 


HD 188603 


11.4 


45.5 


168.6 


2.47 


1000x4 ms 




2x5000x12 ms 


04-11-2004 


02:47 


- 03:54 


HD 25604 


5.1 


60.7 


117.1 


1.75 


1000x4 ms 


16000x12 ms 


2x5000x12 ms 


04-11-2004 


03:54 


- 04:36 


TTauN 




61.6 


114.0 


1.56 


1000x4 ms 


16000x12 ms 




04-11-2004 


04:36 


- 05:04 


TTauS 




62.3 


111.4 


1.47 


1000x4 ms 


16000x12 ms 


2x5000x12 ms 


04-11-2004 


05:04 


- 05:29 


HD 20644 


14.7 


59.0 


101.5 


1.69 


1000x4 ms 


12000x12 ms 


2x3000x12 ms 


04-11-2004 


07:18 


- 07:36 


HD 37160 


6.5 


61.0 


107.4 


1.21 


1000x4 ms 


8000x12 ms 


2x2000x12 ms 


04-11-2004 


09:00 


- 09:23 


HD 50778 


17.3 


61.0 


112.6 


1.04 


1000x4 ms 


8000x12 ms 


2x2000x12 ms 



* determined for the fringe tracking sequence 
** determined for the photometric measurements 
"* peculiar instrumental visibility or chopping problems 



scans are repeated typically 100-200 times in saw-tooth manner 
to increase the statistical accuracy. For typical exposure times 
of 15 ms and readout times of 3 ms this means that the temporal 
fringe signal is modulated at a frequency of « 10 Hz. For these 
inteferometric measurements no chopping is used, because the 
two resulting output signals are subtracted from each other as 
first step of the data reduction. This leads - when supplemented 
with high-pass filtering - to an efficient background subtraction. 

To derive the (wavelength dependent) visibilities, spectra are 
measured with the same pixels of the detector that have been 
used before for measuring the fringe signal. To accomplish this, 
first the light from one, then from the other telescope is blocked. 
For the photometric measurements the secondary mirrors of the 
UTs are chopped to remove the thermal background radiation. 

By definition the visibility is obtained as the ratio of cor- 
related flux (from the interferometric measurement) to the total 
flux (from the spectrophotometry). This gives the raw visibil- 
ity, still suffering from atmospheric and instrumental correlation 
losses. To correct for these losses, the sequence of observations 
has to be repeated for a source with known diameter, and there- 
fore known visibility. 

2.2. Calibration 

Calibrator stars with known diameter were observed immedi- 
ately after the object and in the same region of the sky to correct 
for the reduction in the fringe contrast due to optical imperfec- 
tions and atmospheric turbulence. However, with the present ac- 
curacy of about 5-10% per single visibility measurement, it is 
also possible to use calibrators observed in the same mode dur- 
ing the same night. The calibrators given in Table 1 with their 
HD number were taken from the "MIDI Calibrator Catalogue" 1 
of 509 stars with a flux of at least 5Jy at 10 fim and selected 
for the absence of circumstellar emission, close companions that 



could disturb the visibility measurement, or strong variability. A 
subsample of these calibrators can also be used for absolute flux 
calibration 2 . 

2.3. Additional Photometric Datasets 

In addition to the N-band spectra obtained with MIDI 
(Section 4.2), we derived mid-infrared fluxes from measure- 
ments with the "Infrared Space Observatory" (ISO). For the re- 
duction of the raw data the "ISOPHOT Interactive Analysis" 
software (PIA) version 10.0.0 was used. The long wavelength 
observations were small rasters, producing at 60 and 100//m fi- 
nal maps of 9 x 3 pixels with a pixel scale of 43" x 43" and 
at 150 and 170/im maps of 6 x 2 pixels with a pixel scale of 
89" x 89". The flux distribution of these small maps was fitted 
by a sum of a theoretical PSF and a small extended clump (as- 
sumed in Gaussian shape). This approach was chosen to elim- 
inate the contributions from the nebulae that are present in the 
T Tau system. Details of the method are described in Abraham 
et al. (2000). Nonetheless, due to the large beam size of the tele- 
scope the colour-corrected ISO fluxes listed in Table 2 will be 
considered as upper limits only. 

In the visual, we took images in four bands by using MPIA's 
KING telescope on the Konigstuhl in Heidelberg with its 0.7 m 
primary mirror. These images have been reduced and analysed 
with standard IRAF routines. For the near-infrared, the sub- 
mm, and the mm-regime, fluxes from the literature have been 
adopted, chosing as a rule data corresponding to the low state of 
T Tau Sa. For a discussion of the Q-band fluxes see Section 6.2. 
All photometric measurements are listed in Table 2. 



1 http://www.eso.org/~arichich/download/vlticalibs-ws/ 

2 http://www.eso.org/sci/facilities/paranal/instruments/midi/tools/ 
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Table 2. Photometric Measurements 



Band 


^central 


Date 


Instrument 


TTauN 


T Tau Sa 


T Tau Sb 


Reference 




[/mi] 














TT 

U 


U.JU 




PTir\t / \/l 1 A/TtiirliiniiV 
rllUl. / IVll. IVlalUallaJx 


11 74. + fl 04 mad 






a 




0.43 


Oct 7004 


KTNG / MPT A 


11 7+01 macr 






b 


V 


0.55 


Oct. 2004 


KING / MPIA 


9.9 + 0.1 mag 


_ 


_ 


b 


R 


0.70 


Oct. 2004 


KING / MPIA 


8.9 + 0.1 mag 


- 


- 


b 


I 


0.90 


Oct. 2004 


KING / MPIA 


8.3 + 0.1 mag 


- 


- 


b 


J 


1.25 


14-12-2002 


NACO / VLT 


7.1 +0.1 mag 


> 14.5 mag 


12.9 ± 0.1 mag 


c 


H 


1.65 


14-12-2002 


NACO / VLT 


6.2 ±0.1 mag 


12.3 ± 0.1 mag 


10.4 ± 0. 1 mag 


c 


Ks 


2.16 


14-12-2002 


NACO / VLT 


5.7 + 0.1 mag 


9.2 ±0.1 mag 


8.8 ±0.1 mag 


c 


L' 


3.80 


13-12-2002 


NIRC-2 / Keck 


4.32 + 0.05 mag 


5.64 ± 0.03 mag 


6.25 ± 0.03 mag 


tie 


Ms 


4.67 


13-12-2002 


NIRC-2 / Keck 


4.13 + O.lOmag 


4.96 ± 0.06 mag 


5.77 ± 0.06 mag 


df 




8.5 


Oct./Nov. 2004 


MIDI/VLTI 


5.6 + 0.2 Jy 


2.6±0.2Jy 


1.0±0.1Jy 


b 




o r\ 
y.yj 


Uct./JNOV. ZUU4 


MIDI/ V LI 1 


o. / + v.J Jy 


1 C i A 1 Tit 

l.j ± u.z Jy 


u.o ± u.i Jy 


b 




9.5 


Oct /Nov 2004 


MIDI/VLTI 


7.4 + 0.5 Jy 


1.0±0.2Jy 


0.4±0.1Jy 


b 




10.0 


Oct /Nov 2004 


MID I /VI TT 


8 4 + 4 Tv 


1.1 ±0.2 Jy 


0.5±0.1Jy 


b 




10.5 


Oct./Nov. 2004 


MIDI/VLTI 


8 8 + 4 Jv 


1.6±0.3Jy 


0.7±0.2Jy 


b 




11.0 


Oct./Nov. 2004 


MIDI/VLTI 


8 8 + 4 Jv 


2.3±0.3Jy 


0.9±0.2Jy 


b 




11.5 


Oct./Nov. 2004 


MIDI/VLTI 


8.7 + 0.4 Jy 


3.2±0.4Jy 


1.1 ±0.2 Jy 


b 




12.0 


Oct./Nov. 2004 


MIDI/VLTI 


8.4±0.5Jy 


4.4±0.5Jy 


1.3±0.2Jy 


b 




12.5 


Oct./Nov. 2004 


MIDI/VLTI 


7.8 + 0.5 Jy 


5.3±0.6Jy 


1.3±0.3Jy 


b 


O 

V 


19.91 


27-08-1996 


MAX / UKIRT 


20 + 5 Jy 


28±7Jy 


S 




25 


25-09-1997 


ISOPHOT 




40.88 ± 6.13 Jy 




b 




60 


25-09-1997 


ISOPHOT 




117.77 ± 17.67 Jy 




b 




100 


25-09-1997 

— - ' V / y i y y l 


ISOPHOT 




105.55 ± 15.83 Jy 




b 




150 


25-09-1997 

— ~ < v./ y i y y i 


ISOPHOT 




1 13.39 ± 17.01 Jy 




b 




170 


25-09-1997 


ISOPHOT 




82.79 ± 12.42 Jy 




b 




350 


2004-2005 


CSO 




8. 149 ±0.253 Jy- 




h 




450 


1997-2002 


SCUBA / JCMT 




1.655 ±0.218 Jy 




h 




840 


28-10-1994 


JCMT & CSO 




1.350 ±0.675 Jy - 




i 




850 


1997-2002 


SCUBA /JCMT 




0.628 ± 0.017 Jy- 




h 




1100 


01-12-1994 


OVRO 


0.397 ± 0.035 Jy 


< O.lOOJy 


i 




1300 


April 1988 


IRAM 




0.280 ± 0.009 Jy- 




h 




2800 


1996-1997 


BIMA 


0.050 ± 0.006 Jy 


< 0.009 Jy 


j 



" Melnikov, S. (private communication) 


f Ghez et al. (1991): flux of T Tau N 


* this work 


s Herbst et al. (1997) 


c Herbst et al. (2007) 


* Andrews & Williams (2005) and references therein 


d Duchene et al. (2005): flux ratios 


' Hogerheijde et al. (1997) 


" Beck et al. (2004): flux of 1 lau N 


' Akeson et al. (1998) 



The derived spectral energy distributions (SEDs) are dis- 
played in Figures 4, 12, and 14. The ISO fluxes therein are dis- 
tributed to T Tau N and T Tau S according to the flux ratio found 
in the upper N-band between 1 2 and 1 3 fim. For the flux ratio 
of T Tau Sb and T Tau Sa we assume a value of 0.5. The latter is 
also applied to the Q-band flux of T Tau S. 



3. Data Reduction 

3.1. Visiblities 

For the data reduction, a custom software called MIA+EWS 3 ' 4 , 
written in IDL and the C language, was used, where we chose the 
first branch based on power spectrum analysis (called MIA). This 
package works best when OPD scans are centred on the white- 
light fringe (zero OPD), as in most of our observations, but can 
also handle data tracked at non-zero OPD. The data reduction 
steps are described in a tutorial on the cited web page, as well 
as in Leinert et al. (2004) and Ratzka (2005). The results are 



raw values for the 8-13 /J.m spectrum F mw (A) and the 8-13/Lan 
correlated flux (A). 

To correct for the bias introduced by signal fluctuations not 
related to the fringe signal, the off-fringe power spectrum, deter- 
mined on the source, but far from zero OPD, is subtracted from 
the signal. The fringe amplitude (or correlated flux) as a function 
of wavelength is the square root of the fringe power spectrum 
after this "background" subtraction. Division of this final fringe 
amplitude by the photometric flux gives the raw visibility of the 
object as a function of wavelength. To allow at least approxi- 
mately for the influence of unequal fluxes in the two interfering 
beams, the raw visibilities V raw are calculated by 



Vraw W = 



rcor 
raw 



V*£wU)f£wU)' 



(i) 



where A and B are the fluxes determined from the photometric 
datasets of the two incoming telescope beams. A precise cor- 



3 http://www.strw.leidenuniv.nl/~nevec/MIDI/ 

4 http://www.mpia-hd.mpg.de/MIDISOFT/ 
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Fig. 1. The calibrated visibilities for T Tau N (upper row) and T Tau S (lower row). The errors are the calibration errors. The thick lines within the 
concentric rings on the right side of the individual panels represent the length and position angle of the projected baseline. Between 9.4/jm and 
9.9 /jm the calibrated visibilities for T Tau S are not used for analysis. This wavelength range is strongly affected by the atmospheric ozone band 
due to the low flux of the source in the silicate band. 



rection, which would imply determining the fluxes in the two 
beams for each individual interferometric scan or even each ex- 
posed frame, is not possible with the instrument mode used. 

Calibrated visibilities for an object are obtained by dividing 
its raw visibility by the instrumental visibility derived from one 
or more calibrators with known diameter observed within the 
same night: 



(2) 



The second branch of the data reduction software, called 
EWS, uses a clever shift-and-add algorithm in the complex plane, 
averaging suitably modified individual exposures to obtain the 
complex visibility A(/l)-e"* (/i) (see Jaffe (2004) and the documen- 
tation on the web page just referenced). Thus, this algorithm is 
"coherent". It works best if the OPD scan of the interferomet- 
ric measurements is offset from zero OPD. For the low spec- 
tral resolution of the prism (A/AA « 30) an offset of about 5 A 
is suitable. However, EWS can also be used safely without an 
offset from zero OPD when the observed sources are reasonably 
bright so that the background in no major source of error. We 
thus used EWS to confirm the visibility results obtained with 
the MIA power spectrum analysis, but in particular analysed the 
phases computed by EWS to identify the brighter source in the 



close binary T Tau S (see Section 6.2 and Appendices B & C for 
details). 

3.2. Photometry 

While the determination of the normalised and calibrated visi- 
bilities requires the special software package, the absolute spec- 
tra F (A) can be derived from the photometric measurements 
^raw W an d ^raw W by following the standard procedures ap- 
plied when reducing mid-infrared spectra. To do this, however, 
at least a few of the interferometric calibrators have to also 
be spectrophotometric standard stars with well-known spectra. 
Furthermore, the different airmasses of the scientific targets and 
the calibrators as well as the temporal variations in seeing con- 
ditions have to be taken more seriously than for the visibility 
determination. 

3.3. Correlated Flux 

The wavelength-dependent correlated flux F COT (A) can now be 
derived by a simple multiplication of the above determined 
quantities: 



F COI (A) = V(A)-F(A) 



(3) 
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Although one may think of directly calibrating (A) instead 
of using Equation (3), we determine in this paper the correlated 
flux of our comparatively bright sources with the already cali- 
brated quantities V (A) and F (A). 

Since the correlated spectra reflect to a first appromixation 
the emission from the region not resolved by the interferomet- 
ric instrument, they carry valuable information about the depen- 
dence of the spectrum on the distance from the star. 



4. Direct Results 

In the following, we present the visibilities and photometry of 
T Tau N and T Tau S. These results can be directly derived from 
the observations with MIDI without making assumptions or us- 
ing simulations. 



4.1. Visibilities 



20 



— 10 




10 11 
Wavelength [fj.rn] 

Fig. 2. The spectra for the northern and the southern component to- 
gether with a combined spectrum of the whole source. A spectrum ob- 
tained with Spitzer in February 2004 is drawn as dashed black line. 
Underlaid in grey are two spectra taken with TIMMI 2 in February 2002 
(solid line) and December 2002 (dashed line). 



Separate interferometric observations were possible for T Tau N 
and the T Tau S subsystem. The results are shown in Figure 1 . 
Here, the calibration is performed with all calibrators of one 
night that show a reasonable visibility. Calibrator measurements 
that have not been used are indicated in Table 1 . The errors rep- 
resent the standard deviation of the ensemble of calibrator mea- 
surements. For the two measurements of the scientific targets 
with corrupted photometry (see Table 1) we used the photometry 
for the object in the two other nights and averaged the resulting 
visibilities. In this case, the error bars of the average visibility 
span the whole range that was originally covered by the error 
bars of the visibilities averaged, i.e., the new error bars form an 
"envelope" around the original visibilty errors. 

Both T Tau N and T Tau S are well resolved with visibility 
values corresponding to a FWHM of a Gaussian brightness dis- 
tribution of ~ 10 ± 5 mas at 8 yum. The FWHM roughly increases 
with wavelength to ~ 20 + 5 mas at 1 3 /im. However, the unequal 
appearance of the visibilities of T Tau N and T Tau S calls for 
a different treatment of the two measurements (Sections 5&6). 
While the measurements of T Tau N show a wavelength depen- 
dence of the visibility typical for emission from a circumstellar 
disc (Leinert et al. 2004), the visibilities for T Tau S show sinu- 
soidal oscillations. These have to be interpreted as the signature 
of the binary components T Tau Sa and T Tau Sb. The deep min- 
ima in the modulations already indicate that the brightness of the 
two sources is similar in the N-band. 



4.2. Photometry 

We obtained separate spectrophotometric observations for 
T Tau N and the T Tau S subsystem (Figure 2). From the com- 
parison with TIMMI 2 spectra taken by Przygodda (2004) and a 
Spitzer spectrum, we tentatively conclude that no large bright- 
ness changes occured in the T Tau system between 2002 and 
our measurements. This is similar to what was observed in the 
K-band (Duchene et al. 2005), where T Tau S appears to be in 
a low brightness state. However, it is worth mentioning that the 
spectra deviate significantly from each other, i.e., we witness a 
continuous change of the contributions of the southern compo- 
nents to the overall spectrum. 



5. The northern component 

5.1. Radiative Transfer Model 

We used the radiative transfer code MC3D (Wolf et al. 1999, 
2003; Schegerer et al. 2008) to model simultaneously the spec- 
tral energy distribution and the visibilities of T Tau N by manu- 
ally adjusting the parameters. Manual adjustments have the ad- 
vantage to allow even in the paramater space of our radiative 
transfer model with its many dimensions the determination of 
the best fitting model. On the other hand, this approach always 
leaves a small chance that only a local instead of a global min- 
imum has been found. But we carefully checked this possibilty 
and are confident that we determined the best parameters. 

The used MC3D code includes the effects of heating of a 
disc and an envelope by the central star as well as heating by 
accretion. The photospheric spectrum of T Tau N is represented 
by a Kurucz model (Kurucz 1992) of a star with temperature 

= 5250 K, surface gravity logg = 3.7 [log cm/s 2 ], and solar 
metallicity (Table 3). 

Density Model The density profile of our radial-symmetric 
passively-heated (dust) disc is described in cylindrical coordi- 
nates (r, z) that refer to the distance from the star in the midplane 
and the height above it as 



Pdisc 0,z) = Po(-y) 



3(M) 





z 


H 


[h(r)\ 



)■ 



(4) 



where R+ is the stellar radius and h (r) the disc scale height (Wolf 
et al. 2003). The latter is defined as 



h 0) = h m 



U00AU/ 



(5) 



with /zioo being the scale height of the disc at r = 100 AU. Since 
the Equations (4) and (5) fully describe the structure of the disc, 
po can be adjusted to reach the assumed disc mass. 

As a representation of a potential "envelope", we add a sim- 
ple spherical dust configuration to the disc model. This spherical 
envelope allows us to bring additional material close to the star to 
reprocess light more efficiently. It may represent remnant mate- 
rial of the star formation process. The envelope is geometrically 
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constrained by the inner radius R m and the outer radius R oM of 
the disc. With the density distribution of the disc pd; sc (r) and the 
position vector r, the density profile of the envelope is given by 

/ \r\ \ C2 

PenvO) = Cl • Pdisc(^in, 0) • I — I for R in < \r\ < R out , (6) 

where c\ <K 1 and < 0. While the constant determines 
the concentration of the dust towards the star, the quantity c\ is 
adjusted to guarantee a low optical depth of the envelope and the 
possibility to observe the innermost region of the disc. The disc 
and envelope are combined by 

„ ( r \-\ Pdisc W for p env (r) < PdiscO) 

PenvW - | pem(r) for pmv(r) > pd . sc(r) , V) 

ensuring a smooth transition from the disc to the envelope. 

Accretion Effects We extend the passive model further for the 
implementation of accretion effects (Schegerer et al. 2008). 
Apart from the parameters of the disc and the star, this accretion 
model requires three additional parameters: the accretion rate M, 
the radius /?t,nd at which the disc is magnetically truncated, and 
the boundary temperature 7t, n d of the accreting regions on the 
surface of the star. Following the approach of Konigl (1991) we 
assume here that the truncation radius is half the corotation ra- 
dius. A stellar magnetic field of only about 0.4 kG is then re- 
quired (Johns-Krull 2007). For such a steady-state scenario and 
under the assumption that the shock front is located close to 
the stellar surface, 7t, n d can be easily calculated when the stel- 
lar mass, the stellar radius, the stellar angular velocity, and the 
accretion rate are known (Konigl 1991). With a rotation period 
of 2.8 days (Herbst et al. 1986), we find for the boundary temper- 
ature 7b n d of T Tau N the value listed in Table 3. Interestingly, 
the magnetic truncation occurs at a radius of /? Dn d = 5.4 R in 
this accretion model. The boundary radius is thus smaller than 
Ri n = 23 R Q and the magnetic truncation has to take place in a 
gaseous inner disc. 



a dust mixture of "astronomical silicate" and graphite with rel- 
ative abundances of 62.5% for astronomical silicate and 37.5% 
for graphite (Draine & Malhotra 1993). We use the optical data 
of "smoothed astronomical silicate" and graphite published by 
Draine & Lee (1984) and Weingartner & Draine (2001) and 
consider a grain size power law n(a) oc gt 3 - 5 with a m [„ ^ a ^ 
am, where n(a) is the number of dust particles with radius 
a. This power law was found by Mathis et al. (1977), here- 
after MRN, in their study of the extinction of interstellar carbon 
and silicate with typical sizes between 0.005/im - 0.01/im and 
0.025 yum - 0.25/im, respectively. This grain size power law has 
already been used in former modelling aproaches of discs around 
young stars, e.g., Wolf et al. (2003). We use a minimum particle 
size of a m ; n = 0.005 yum. 

Furthermore, we implement a two-layer dust model. The disc 
interior contains a maximum dust grain size of 1 mm while the 
MRN grain size distribution with a max = 0.25 fim is used in the 
disc atmosphere. These upper disc layers are defined by an op- 
tical depth tn in the N-band, measured vertical to the disc mid- 
plane, below unity. Such a division of the disc is based on the 
idea of the favoured settling of larger dust grains, e.g., Schrapler 
& Henning (2004). Similar models with different grain sizes 
have been investigated for other discs by D' Alessio et al. (1997), 
Men'shchikov et al. (1999), and Wolf et al. (2003). 

In order to avoid determining the temperature distribution of 
each single dust component and to accelerate the radiative trans- 
fer simulations, we average the optical properties of carbon and 
astronomical silicate for different grain sizes in each dust layer 
to derive the optical constants. Such an approach was justified 
by Wolf (2003). 

Interstellar Extinction Since the interstellar extinction is not 
taken into account by the model itself, we "redden" the fluxes 
resulting from the simulation before comparing them with the 

Table 3. The parameters of the radiative transfer model of T Tau N. 



Dust Model Following Schegerer et al. (2008), we assume in 
our modelling approach a gas-to-dust mass ratio of 100:1 and 



Parameter 



Unit 



Value 



Ref. 




0.1 1.0 10.0 100.0 1000.0 

Wavelength [jim] 

Fig. 3. Comparison of the mass absorption coefficients of the dust in 
the disc atmopshere (solid) with those of the insterstellar extinction law 
(diamonds) from Rieke & Lebofsky (1985). The values are given rel- 
ative to the V-band. Plotted as dotted curve are the opacities provided 
by Ossenkopf & Henning (1994) for a gas density of 10 6 cirr 3 , thin ice 
mantles, and a coagulation over 10 5 years. 



Stellar Parameters 

Mass (M*) M 

Temperature (7\) K 

Luminosity (L*) L Q 

Radius (i?,) R Q 

Accretion 

Accretion Rate (M) M Q yr~ 

Boundary Temp. (T bni ) K 

Boundary Rad. (Rbnd) R* 

Circumstellar Disc / Envelope 

Inner Radius (R in ) AU 

Outer Radius (R out ) AU 
P 

h m AU 

Cl 

Inclination deg 

Disc Mass (M disc ) M Q 

Interstellar Extinction 

Foreground Extinction A v mag 



2.1 
5250 
7.3 
3.3 



•10" 8 
9700 
1.6 



0.1 
80 
1.25 
18 

■ io- 5 

-5.0 

< 30 

■ io- 2 



1.5 



calc. 



calc. 
calc. 



a White & Ghez (2001); * this work; c Hogerheijde et al. (1997); 
calc. = calculated 
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Fig. 4. Results of the simultaneous fit to the SED and the visibilities of T Tau N. Top Left: The model (Table 3) is plotted for inclinations of 
10°, 20°, and 30°. The contributions from the stellar photosphere based on a Kurucz model (dashed) and the accretion (dotted) take foreground 
extinction into account. Measured photometric data are indicated by black crosses. Arrows indicate upper limits. For references, see Table 2. Top 
Right and Lower Row: The measured visibilities (error bars) are displayed together with the upper and lower limits (squares and triangles) derived 
from the model with an inclination of 20° by varying its orientation by 360°. 




-100 -50 50 100 -100 -50 50 100 -100 -50 50 100 

u [arcsec -1 ] u [arcsec -1 ] u [arcsec -1 ] 



Fig. 5. The modelled visibilities of T Tau N for 8.4 pm, 10.6 pm, and 12.6 pm and a disc inclination of 20° (from left to right). The solid contours 
are drawn for steps of 0. 1 and the associated grey background is underlaid for easier comparison. The radii of the overplotted dashed circles 
correspond to the spatial frequencies for which the baselines are sensitive at the various wavelengths. 



measured values. Starting with the visual extinction Ay the ex- 
tinctions for the other wavelengths are calculated by applying 
the wavelength-dependent mass absorption coefficients. These 
or more presisely the cross sections are provided by the model of 
the dust in the disc atmosphere described in the last paragraph. 
A comparison with the interstellar extinction law by Rieke & 
Lebofsky (1985) is shown in Figure 3. The visual extinctions Ay 



of the models are given together with the other parameters in 
Tables 3 & 5. 

5.2. Results of Model Fits 

The parameters of our best model for T Tau N are listed in 
Table 3. In Figure 4, the resulting SEDs and normalised visibil- 
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1000.0 



Fig. 6. The spectral energy distribution for the model of T Tau N with a 
disc inclination of 20° when neglecting the accretion (solid), removing 
the "envelope" (dotted), and changing the outer radius of the disc to 
50 AU (dashed), 110 AU (dashed-dotted), respectively. 



ities are plotted together with the measured data. A comparison 
shows that the model and the data are in good agreement, al- 
though the stellar parameters, the disc mass 5 , and the accretion 
rate were taken as constrained by complementary studies. 

However, there are differences between the measured SED 
and the model. The slight deviations in the N-band are con- 
sequences of the simplistic dust model. The "underestimation" 
of the measured SED in the wavelength range from 25 fim to 
170/im can be attributed to the large beam size of the ISO satel- 
lite. We already treated the measured fluxes as upper limits, be- 
cause the nebulae around T Tau N contribute significantly at 
those wavelengths. The same might be applicable for the Q- 
band. Here, the resolution of the telescope was limited and the 
southern binary could only marginally be resolved (Herbst et al. 
1997). 

Due to the uncertainties at wavelengths longwards of the N- 
band, models with various combinations of the parameters /? and 
/iioo fit the SED and the visibilities. Also R out is not well con- 
strained. Figure 6 shows the changes in the SED when vary- 
ing the outer radius of the dies R out by 30 AU while keeping 
all other model parameters constant. The visibilities both in the 
near- and the mid-infrared are only marginally affected. In a re- 
cent study, based on the extent of the H2 emission and models 
of the SED, an outer radius of the disc of about 85-100 AU was 
found (Gustafsson et al. 2008). The value listed in Table 3 thus 
probably underestimates the true size of the disc. 

Nonetheless, our model puts strong constraints on the incli- 
nation of the disc. Large inclination angles can be safely ex- 
cluded on the basis of both the SED and the visibilities. The 
different levels of the measured visibilities only reflect the spa- 
tial resolution reached by the different baseline configurations 
and are not an indication for an elongated disc (Figures 4&5). 
On the other hand, since the variations of the visibilities due to 
the variation of the position angle are small, the determination 
of the actual position angle of the major axis of the disc around 
T Tau N is not possible. 

Another constraint of the model is the presence of a thin 
envelope with a steep density gradient, where the material is 



highly concentrated towards the central star. Without this adi- 
tional material the flux between 3 /im and 8 fim would be too 
low (Figure 6) and the visibilities represent no longer the mea- 
surements. The visibilities are reduced by about 50% when com- 
pared to our best model. However, it has to be mentioned that 
models with a curved inner rim or other special geometries of the 
inner disc have not been tested, although such models might not 
require the steep density gradient of the envelope found above. 

Although accretion may lead to similar effects as an ex- 
tended, spherical dust distribution (Schegerer et al. 2008), ac- 
cording to the measurements the accretion rate of T Tau N is 
not high enough. However, our model is not very sensitive to 
the accretion rate. Neither the visibilities nor the SED change 
significantly (Figure 6). 

5.3. K-band visibilities 

The interferometric measurements in the mid-infrared presented 
in this paper constrain mainly the structure of the warm inner 
parts of a disc, whereas for the study of the hot innermost regions 
near-infrared visibilities are required. T Tau N was measured 
with the "Palomar Testbed Interferometer" (PTI) at 2.2 fim both 
with the north-south (103-108m) and the east-south (83-86m) 
baseline (Akeson et al. 2002). After correcting for the incoherent 
contribution of T Tau S (Akeson et al. 2000) the squared visibil- 
ities are roughly between 0.6 and 0.7 for the short and 0.4 and 
0.5 for the long baseline. For our best-fit model (Section 5.2) 
we find at 2.2 fim squared visiblities of 0.59 + 0.01 for a 85 m 
baseline and 0.45 + 0.01 for a 105 m baseline. The errors given 
here reflect the uncertainty due to the unknown position angle 
of the disc. These predictions of our model are in good agree- 
ment with the PTI measurements, although no special attempt 
has been made to fit the near-infrared results. 



5.4. Dust Distribution and Evolution 

For the observations on the VLTI, T Tau N is well separated from 
T Tau S in the 10 /im range, so we can study its uncontaminated 
mid-infrared spectrum. The modest extinction (Table 3) and the 
fact that the silicate band appears in emission allows us to in- 
terpret the profiles by a mixture of amorphous and crystalline 
particles of various sizes. By applying such studies to the total 
and correlated spectra of selected Herbig Ae/Be stars, evidence 
for grain growth and crystallisation particularly in the inner parts 
of the circumstellar discs was shown (van Boekel et al. 2004). 
Similar processes in the circumstellar discs of T Tauri stars were 
expected and recently found for RY Tau (Schegerer et al. 2008) 
and even for the less luminous source TW Hya (Ratzka et al. 
2007). 

Here, we want to analyse the total and correlated spectra of 
T Tau N (Figure 7) to determine the silicate composition of its 
circumstellar disc down to radii of a few AU. For this we use 
a ^-fitting method (Bouwman et al. 2001). It assumes that the 
silicate emission feature has its origin in the optically thin sur- 
face layer of the circumstellar disc, where it results from a linear 
combination of mass absorption coefficients of different dust 
species ;' of different sizes j: 



F(v) = B(v, T) 



5 The opacities used by Hogerheijde et al. (1997) are plotted for com- 
parison in Figure 3 as dotted curve (M. Hogerheijde, priv. communica- 
tion). Since MC3D only provides cross sections no absolute opacities 
can be given. 



Co + J] C U "ijW 



(8) 



where Co and C,y are fitting parameters, which reflect the under- 
lying continuum and the mass contribution of each component. 
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Fig. 7. Contributions of the dust species to the total flux of T Tau N (upper left panel) and the correlated fluxes. The three different baselines 
are indicated. In each panel the best fit (dashed black) is plotted together with the measured spectrum (solid black). The lines above the fitted 
continuum (solid grey) represent the various contributions: small amorphous (dotted grey), large amorphous (dashed grey), and crystalline grains 
(dashed-dotted grey). 



The quantity F(v) is the spectral flux at frequency v, /Q/(v) rep- 
resents the frequency-dependent mass absorption coefficient for 
a specific component, and B V (T) is the Planck function corre- 
sponding to a blackbody temperature T. 

Observationally, the signature of pristine, interstellar dust 
particles is a broad bell-shaped emission feature with a maxi- 
mum near 9.7 /im. Larger particles are characterised by a flat- 
topped emission feature, while crystalline particles show nar- 
rower emission, most notably one at 11.3yum due to forsterite. 
As basic dust set for our ^-fitting routine we use the same 
silicate species already used by Schegerer et al. (2006): small 
(0.1 fim) and large (1.5 fim) amorphous grains with olivine and 
pyroxene stoichiometry (Dorschner et al. 1995), as well as crys- 
talline species such as forsterite (Servoin & Piriou 1973), en- 
statite (Jager et al. 1998), and quartz (Henning & Mutschke 
1997). Carbon is not considered in the study presented here, be- 
cause its emission profile is monotonic in the 10 /im wavelength 
range (Draine & Lee 1984; Wolf & Hillenbrand 2003) and thus 
contributes to the underlying continuum only. 

In the case of T Tau N, we find from fitting the emission pro- 
file of the total flux, corresponding to an emitting region of about 
1-20 AU, that small and large amorphous grains dominate with 
an abundance of about 90% (Table 4). On the other hand, the 
correlated flux spectra clearly show the dominant role of large 



amorphous and crystalline species (see Figure 7). These inter- 
ferometric measurements allow insights into the dust properties 
at radii close to the central source. The nominal resolutions A/B 
of the three baselines are 7AU, 5 AU, and 4AU at 10/im and 
the distance of T Tau. Here, the abundance of small amorphous 
grains is insignificant. 

However, the trend for an increasing contribution from crys- 
talline species with spatial resolution is broken in the case of the 
62 m baseline. Here, one has to keep in mind that the interfer- 
ometer resolves the disc only in the direction of the projected 
baseline, while perpendicular to it, light is also collected from 
outer disc regions. The determination of the absolute positions 
and abundances of the dust species thus requires a more detailed 
study. Nevertheless, the finding that large amorphous and crys- 
talline grains are concentrated towards the central source appears 
to be robust. 



6. The close binary T Tau S 

The observations of T Tau S are complicated by the fact that its 
two components are not resolved in the N-band by a single VLT 
telescope, and are therefore superposed both in the photometric 
and in the interferometric measurements. In the latter the effect 
of binarity is visible as sinusoidal modulations in the visibilities 
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Table 4. Results of the spectral decomposition of the MIDI spectra. 
Given are the mass fractions of the various components (rows 1, 2, and 
3) and the temperature of the continuum (row 4) for the total spectrum 
and the correlated spectra. The first two rows give the values for the 
small and large amorphous dust gains. The third row lists the contri- 
bution from both the small and large crystalline grains, i.e., enstatite, 
forsterite, and quartz. 



Species 


total 


43 m 


62 m 


85 m 


0.1 yum 


0.24 ±0.14 


0.00 ± 0.00 


0.00 ± 0.02 


0.00 ± 0.00 


1.5 ^m 


0.64 ±0.10 


0.63 ± 0.05 


0.80 ± 0.03 


0.59 ± 0.24 


cryst . 


0.12 ±0.03 


0.37 ± 0.07 


0.20 ± 0.02 


0.41 ± 0.06 


Temp. 


355 ±5K 


403 ± 14 K 


453 ± 4 K 


490 ± 25 K 



(Figure 1). In the following, inspired by the procedure routinely 
used in speckle interferometry, we try a decomposition in order 
to be able to study the two components individually. 

This requires several steps. First of all, we determine from 
the interferometric measurements the binary parameters, i.e., 
separation, position angle, and the ratio of the correlated fluxes 
(Section 6.1). Then the brighter component is identified by the 
measured phases. After estimating the visibilities of the fainter 
component by means of a reasonable radiative transfer model of 
its SED, we use these visibilities to calculate separate N-band 
spectra as well as visibilities for T Tau Sa (Section 6.2). Finally, 
a radiative transfer model for T Tau Sa is presented (Section 6.3). 

6.1. Binary Parameters of T Tau S 

All visibilities of T Tau S show that this source is well resolved 
by the interferometer (Figure 1). In addition, a clear sinusoidal 
binary signal is apparent for the second and third night when the 
position angle of the projected baseline was fairly close to the 
position angle of the separation vector of the two components. 
In Figure 8 the calibrated visibilities are plotted as functions of 
the spatial frequency u = B/Afor ease of discussion (compare 
Figure 1). Here, B is the length of the projected baseline. 

In Appendix A, we derive a formula for the visibility of a 
binary composed of two extended components. Under the as- 
sumption that the difference of the Fourier phases is negligible 
one finds: 



V(u) 



F c ° r (u) z +2F c ° r (u)F c c ° r (u) cos (2ttus) + F c c 0I (u)' 



**»(«) + Ff\u) 



(9) 



The suffixes p and c refer to the primary and the companion, and 
s is the separation of the two components along the projected 
baseline. We use the definition of the flux ratio 



/(«) 



F™(u) 
Ff(u) 



< 1 



and obtain 



Vfl,(ii) = V («) 



Vl + f 2 (u) + 2f(u) cos (2nus) 
!+/(«) 



(10) 



(11) 



Here, Vq(u) represents the visibility as it would be produced in 
a scan perpendicular to the separation vector, i.e., without the 
binary modulation. The spatial frequency u is here a function of 
the wavelength, while the projected baseline length is fixed. The 
total flux is thus also a function of u. 
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Fig. 8. The calibrated visibilities vs. spatial frequency are given by the 
error bars. Overplotted are the best fits (diamonds) for the measure- 
ments with the two long baselines. The lines are the linear approxi- 
mations of the extended structures surrounding the components of the 
southern binary T Tau S. A fit for the measurement with the short base- 
line has not been attempted due to the limited modulation of the signal. 
The dashed line in the first panel shows the visibility as it is expected 
from the fitted parameters. A factor 0.6 has been applied. 



Fitting the visibility instead of fitting the correlated fluxes 
has the advantage that we can use the simple linear relation 



Vq(u) = do +CL\U 



(12) 



as a reasonable approximation. This is supported by the obser- 
vation of the first night, scanning approximately perpendicular 



12 



Th. Ratzka et al.: Spatially resolved mid-infrared observations of the triple system T Tauri 



to the separation vector. The derived visibility shows a smooth, 
almost linear trend with wavelength or spatial frequency, respec- 
tively. We expect similar smooth trends to underly on our other 
observations. 

The correlated flux ratio of the two components is modelled 
by a second-order polynomial in u as 



/(") = /o + fiu + fiu 2 , f(u) < 1 



(13) 



to allow for the effects of silicate absorption, while the separa- 
tion of the components is formally allowed to vary linearly with 
spatial frequency to take into account that in general the photo- 
centre of an inclined star-disc system can vary with wavelength: 



s = s(u) = sq + S\u , 



(14) 



A non-linear least squares fit, based on the Levenberg- 
Marquardt algorithm, of the visibility function to the visibilities 
measured with the two longer projected baselines has been per- 
formed to derive two independent parameter sets (Figure 8). 6 We 
find rather small values for s\ in both fits, which shows that the 
wavelength-dependence of the separation can be neglected. 7 We 
take this as indication that our ansatz of Equation (1 1) is mean- 
ingful. 

The fit to a visibility described above yields s(u), i.e., the 
separation of the binary projected along the position angle of 
the baseline projected onto the sky. This means that we can only 
derive one component of the two-dimensional separation vector 
from each visibility. Fortunately, we have two fitted measure- 
ments at different orientations of the baseline to determine the 
relative position of the components at the time of the observa- 
tions. 

The two derived projected separations are indicated by the 
dashed lines in Figure 9, where the primary T Tau Sa is located 
at the origin of the coordinate system. The possible positions 
of the companion (black ellipses) are determined by the inter- 
sections of the dashed lines. Two symmetrical solutions exist, 
because Equation (11) does not allow the determination of the 
sign of s(u). Based on observations of the T Tau S binary in the 
near-infrared we chose the western position (solid ellipse). At a 
wavelength of 10.5 jim we thus find 



d = (124 + 8) mas 
= (300 + 5)° 



(15) 



for the separation d and the position angle 0. 

Figure 9 shows that this result fits reasonably well to earlier 
measurements and to the orbit of the T Tau S binary determined 
by Kohler (2008). 

Besides the information on the relative position of the com- 
panion, the sinusoidal modulation of the visibility as a function 
of spatial frequency u (Figure 8) carries information on the flux 
ratio of the components. In detail, the observed visibility mod- 
ulations for T Tau S deviate from a sinusoidal shape, and it 
is this deviation which carries information on the wavelength- 
dependence of the brightness ratio in the 8-13/im range. The 
fit procedure actually takes this into account and predicts the 
brightness ratio as a function of spatial frequency f(u) to have a 



6 Since the separation vector of the binary is almost perpendicular to 
the position angle of the shortest projected baseline, the binary modu- 
lation of the visibility measured with this baseline is not pronounced. A 
fit for the visibility therefore has not been attempted in this case. 

7 Between 8 /mi and 13 /mi, the separation increases from 120 mas 
to 125 mas for the intermediate and from 102 mas to 103 mas for the 
longest baseline. 
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Fig. 9. In grey the near-infrared measurements of the positions of 
T Tau Sb with respect to T Tau Sa are plotted, see Kohler et al. (2008) 
and references therein. The ellipses represent the errors. The best-fit or- 
bit derived by Kohler (2008) is also shown (dotted: periastron, dashed- 
dotted: line of nodes). The positions marked by crosses are separated 
by one year. Overplotted are the projected separations derived from our 
mid-infrared interferometric measurements and the determined relative 
position (black ellipse). The dashed, black ellipse is the second possible 
position due to the 180° -ambiguity. 



maximum (in the sense of being close to 1) near 9.5 /im and min- 
ima of the order of 0.2 at short and long wavelengths (Figure 10). 

Although this separation of the fluxes looks elegant, it is an 
approximation only. What is actually determined and written in 
Equation (1 1) in form of the polynomial f(u) is not the flux ra- 
tio, but rather the ratio of the correlated fluxes. Only with the 
additional assumption that the correlated fluxes of the two com- 
ponents have a similar dependence on the spatial frequency u, 
i.e., the relative light distributions are the same, this ratio also 
represents the ratio of the total fluxes. We use this assumption 
to determine below the separate N-band spectra of T Tau Sa and 
T Tau Sb. 

To put this assumption on a firmer basis, we repeated the 
binary fitting for the correlated fluxes of T Tau S for which 
we describe the "binary-free" correlated flux by a second- 
order polynomial instead of the mere linear pre-factor Vo(u) in 
Equation (12). The second-order polynomial is necessary to fit 
the silicate absorption feature. The relative position found by 
these fits is fully consistent with the position determined above 
from the visibility fits. The ratio of correlated fluxes, indepen- 
dent of the baseline length, is very close to that derived from the 
visibility fits (Figure 10). 

The average of the flux ratios determined from the fits of 
the visibilities and the fits of the correlated fluxes has been 
used to calculate the spectra shown in Figure 1 1 . We consider 
these fluxes as already good estimates, which will be refined in 
Section 6.2. 

Both components of T Tau S have the silicate band in ab- 
sorption, but the optical depth, as usually defined by 



F v = FT 



(16) 
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Fig. 10. The resulting flux ratio of the correlated fluxes when averaging 
the two polynomials derived from the fits to the two visibilities (black). 
Also shown is the flux ratio when averaging the two polynomials de- 
termined by the fits to the two correlated fluxes (grey). The latter are 
shifted by 0.05 fim to longer wavelengths to allow a better comparison 
of the results. 
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Fig. 11. The total spectrum of T Tau S (grey) as well as the individual 
fluxes of T Tau Sa and T Tau Sb under the assumption that the average 
of the flux ratios given in Figure 10 is the flux ratio of the total fluxes. 
Despite this approximation, the results are close to the final spectral 
decomposition (Figure 13). 



is quite different. One finds towards T Tau Sa an optical depth 
of t = 1.8, while towards the Ml component T Tau Sb the value 
is only t = 0.9. For the determination of the optical depth it is 
assumed here that the flux values at 8yum and 13/Lan represent 
the continuum, that this continuum can be linearly interpolated, 
and that the minimum of the silicate absorption is reached at 
9.7 fim (Figure 11). 

Given the small separation between the two components 
(0.1" or 18 AU), the difference in extinction is hard to explain 
by variations in the density of foreground material. We con- 
clude that T Tau Sa probably has additional intrinsic extinc- 
tion in addition to an extinction common to both components. 
Based on near-infrared observations and after a detailed analysis 
Duchene et al. (2005) already found evidence that T Tau Sa is 
surrounded by a small edge-on disc causing the intrinsic extinc- 
tion. If T Tau Sb also possesses a circumstellar disc, it could not 
be aligned to the disc of T Tau Sa, because the extinction then 



would be much too high. Skemer et al. (2008) confirm this mis- 
alignment observationally based on photometric measurements 
in the N-band. Perhaps the circumstellar disc around T Tau Sb is 
coplanar with the orbit of this source around T Tau Sa, which is 
inclined by about 35° (Section 7). 

6.2. Disentangeling the Visibilities of the Two Southern 
Components 

We have the following situation: we measured the visibility 
V m eas (Figure 1) and the total flux F l ^.. ds (Figure 2) of the south- 
ern binary. As explained in Appendix A the normalised visibil- 
ity perpendicular to the separation vector of the binary is the 
flux weighted average of the visibilities of the components, see 
Equation (A.2). This visibility must be close to the slowly de- 
creasing visibility measured for the shortest baselines with val- 
ues between 0.5 and 0.6. By multiplying this "binary-free" visi- 
bility with the measured total flux one can derive the correlated 
flux F™. ds and with the flux ratio / < 1 (Figure 10) the individual 
correlated fluxes F C p r and F™ 1 . To derive the individual visibili- 
ties, the following system of equations has to be solved, wherein 
bold quantities are known 8 : 

pcor _ \r Tjitot 
1 meas " meas ' r me as 



pcor 

meas 



pcor 

meas 



1/(1 +f) 



f/(l+f) = 



pcor 



pcor 



V P 



= V c 



r P 



F tot 
r C 



(17) 



To accomplish this, we will constrain V c , but we first have 
to determine which of the components of T Tau S is the primary, 
i.e., is brighter in the N-band, and whether indeed T Tau Sa is 
the source with the stronger silicate absorption as has been spec- 
ulated above. Intuitively, we would argue that the eastern, the 
more deeply embedded source has to be the brighter component 
in the mid-infrared, but so far we do not know. A phase measure- 
ment is needed. 

Phase values + indicate asymmetry of the object. A typi- 
cal application of phase measurements in a binary is to resolve 
the right-left ambiguity of the position of the companion with 
respect to the primary remaining after analysis of the visibility 
modulus. With the centre of light at the origin, the phase of a 
binary shows a staircase appearance, rising in the direction to- 
wards the companion (if the +x and +u coordinates are going 
in the same direction). If a linear trend in the phase is removed, 
corresponding to shifting the primary to the origin, the binary 
phase will oscillate around zero, rising at the spatial frequency 
of the visibility maxima if a scan over the binary in the +x direc- 
tion would hit the companion first. Our data reduction package 
EWS gives this kind of phase function. We want to use the MIDI 
observations of T Tau S to decide whether it is the eastern com- 
ponent (the infrared companion) or the western component (the 
Ml star) that appears as primary with the larger correlated flux. 

This will only work if the MIDI observing procedure defines 
an orientation on the sky. And indeed, during the OPD scan per- 
formed for interferometric measurements with MIDI, the posi- 
tion corresponding to the white light fringe moves on the sky 



8 Mathematically we cannot exclude that one of the components is 
unresolved. But we expect that this is not the case and that the visibil- 
ity of the faint component is smooth and monotonous, with no obvious 
maximum around 9.5 fim (see the smooth trend of the visibilities ob- 
served for T Tau N and of the visibility measured perpendicular to the 
line connecting the components of T Tau S). 
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Fig. 12. Results of the fit to the SED of T Tau Sb. The model (see text) is 
plotted for inclinations of 30°, 50°, and 70°. The contributions from the 
stellar photosphere (dashed) and the accretion (dotted) take foreground 
extinction into account. Measured photometric data are indicated by 
black crosses. Arrows indicate upper limits. For references, see Table 2. 



along the projection of the baseline in a defined direction: closer 
to that telescope that feeds the right entrance window of MIDI 
from the viewpoint of the incoming light. This direction can be 
determined from the header information given with the data. 

In Appendix B, we derive that the phases, as determined with 
MIDI, indeed should be falling with wavelength (rising with spa- 
tial frequency) at the position of the visibility maxima when the 
companion is offset from the primary component in the direction 
of negative x or OPD, which is the direction defined above. 

In the case of the observations of T Tau S on October 30, 
2004 and November 1, 2004, the phase was rising in the visibil- 
ity maxima, which means that the companion came first when 
moving along the above mentioned reference directions towards 
a position angle of 48° and 85°, respectively, in the sky. The op- 
posite was true on November 3, 2004, when the direction on the 
sky was towards a position angle of 291°. The conclusion from 
each of these observations is that at 10 fim the eastern source, 
the deeper embedded companion T Tau Sa, dominates the corre- 
lated flux and has to be considered as the primary further on in 
this paper. 

The visibilities of the fainter component in the mid-infrared 
Vc are now obtained from a "auxiliary" radiative transfer model 
fit to the measured SED of T Tau Sb. 

Following Duchene et al. (2002) and Duchene et al. (2005), 
we consider T Tau Sb as a T Tauri star of spectral type Ml behind 
an extinction screen of Ay =15 + 3 mag. Using the extinction 
law of Rieke & Lebofsky (1985) and taking the nearby M0. 5 gi- 
ant 6 Oph (HD 14605 1) as a spectral template, we see an infrared 
excess in T Tau Sb increasing by nearly 3 mag from H to L'. The 
near-infrared brightness is probably already affected somewhat 
by additional circumstellar emission. Taking the K brightness of 
T Tau Sb and extrapolating into the mid-infrared with the spec- 
trum of HD 14605 1 gives an upper limit to the photospheric con- 
tribution of T Tau Sb. This upper limit falls from 1 90 mJy at 8 fim 
through 140mJy at 10//m to 85 mJy at 13 /im. These values are a 
factor of 2-10 less than the correlated flux attributed to T Tau Sb 
from our measurements. Obviously, there is strong circumstel- 
lar emission filling up this otherwise unexplained mid-IR flux, 
which we model to come from a circumstellar disc at moderate 
inclination. 



Due to the close periastron passage of the companion, a trun- 
cated disc with a size of not more than 5 AU is expected. Placing 
such a disc with an inner radius of 0.1 AU, an outer radius of 
5.0 AU, /3 = 1.00, /iioo = 7.0 AU, a moderate inclination angle, 
and a mass of 2.5 • 10~ 4 M o around a Ml -type star (Duchene 
et al. 2002), whose photosphere is represented by a black body 
of 3720 K (Schmidt-Kaler 1982), and applying a foreground ex- 
tinction of 15 mag in the visual (Duchene et al. 2005) leads to 
a reasonable representation of the spectral energy distribution of 
T Tau Sb (Figure 12). For this model we adopted a stellar lu- 
minosity of 1.8 L Q (Schmidt-Kaler 1982; Bessell & Brett 1988), 
a stellar mass of 0.4 M Q (Kohler 2008), and an accretion rate 
of 1 • lO^Moyr 1 (Duchene et al. 2005). An "envelope" with 
ci = 5 • 10~ 6 and = -1.0 is also taken into account. It has to 
be mentioned that according to the model the large foreground 
extinction is responsible for the silicate absorption. The disc by 
itself would exhibit the silicate band in emission. 

We refer to Section 5.2 for the explanation of the differ- 
ences between modelled and measured fluxes longward of the 
N-band. But we want to return to the underestimation of the 
Q-band flux. It is not known whether the Q-band flux includes 
contributions from the northern component or the surround- 
ing nebulae, or whether the flux ratio 0.5 between T Tau Sb 
and T Tau Sa as prediced by us based on the N-band fluxes 
is correct. Nonetheless, the main reason for the discrepancy is 
that the southern binary was much brighter in the mid-infrared 
at the time when the Q-band flux was measured. This can be 
seen, when comparing the integrated fluxes of T Tau S reported 
by Herbst et al. (1997) with our results presented in Table 2. 
Depending on the wavelength, the flux was lower by a factor of 
3 to 8 in 2004. 

Due to the uncertainty what fractions of the Q-band and ISO 
fluxes originate in the T Tau Sb system, the model is based on 
the near-infrared and N-band fluxes and is thus ambiguous. In 
the near-infrared the stellar photosphere represents the photo- 
metric measurements always well as long as the inclination of 
the disc is moderate. The N-band fluxes are sensitive to the ac- 
cretion rate, the disc mass, and the inclination angle, but several 
combinations lead to the same quality of the modelled SED. The 
main restriction for a moderate inclination angle comes from the 
near-infrared. Unfortunately, the model does not well represent 
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Fig. 13. The total spectrum of T Tau S (grey) as well as the individual 
fluxes of T Tau Sa and T Tau Sb as determined with the help of the 
modelled visibilities of T Tau Sb. Overplotted as dashed lines are fits 
for astronomical silicates with optical depths of T Sa =l-7 and r S b=0.8. 
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Fig. 14. Results of the simultaneous fit to the SED and the visibilities of T Tau Sa. Top Left: The model (Table 5) is plotted for inclinations of 70°, 
72°, and 74°. The contributions from the stellar photosphere (dashed) and the accretion (dotted) take foreground extinction into account. Measured 
photometric data are indicated by black crosses. Arrows indicate upper limits. For references, see Table 2. Top Right and Lower Row: The black 
squares and triangles represent the upper and lower limits of the normalised visibilities when the orientation of the model for T Tau Sa with a 
disc inclination of 72° is varied by 360° . The grey squares and triangles are the limits for the modelled visibilities of T Tau Sb. The grey dashed 
lines mark the range of values adopted to derive the visibilities of T Tau Sa, which are plotted as solid error bars. The dotted and shifted error bars 
represent the visibilities derived for T Tau Sa, when the visibilities of T Tau Sb are fixed to 1.0, 0.6, respectively, with an error of 5%. The grey 
area between 9.4/jm and 9.9pm marks the atmospheric ozone band that affects the measurements due to the low flux of the sources in the silicate 
band. 



the silicate absorption feature, but a higher foreground extinction 
is excluded again by the near-infrared fluxes. 

However, although the SED is not very well known and the 
model leaves ambiguities, the constraints on size and proper- 
ties of the source suggest that the object visibilities cannot have 
values very different from the derived visibilities (Figure 14). 
This was confirmed by comparing the results obtained for sev- 
eral model assumptions. 

Finally, the upper and lower limits of the modelled visibil- 
ities for a disc inclination of 50° when changing the position 
angle of the disc by 360° have been used to derive the "object" 
visibilities of T Tau Sb. This was done by calculating the aver- 
age and the standard deviation of the limits and interpolating the 
results over the whole N-band (Figure 14). The visibilities found 
for T Ta Sb are typical for a circumstellar disc. 

After determining with the radiative transfer model limits for 
the visibilities of T Tau Sb, we can solve Equations (17) for the 
remaining quantities. To this end we divided the measured visi- 
bilities V me as by the binary signal and multiplied the results with 
the measured spectrum of the southern binary F^ as . The result- 



ing, corrected correlated fluxes of the whole system F™ r as can 
be assigned to the individual components, when taking into ac- 
count that T Tau Sa is the brighter component as has been proven 
by the phase measurements and that the flux ratio / of the corre- 
lated fluxes can be described by the relation plotted in Figure 10. 
Afterwards, we determine the flux of the companion from its 
correlated fluxes F c or with the help of the modelled visibilities 
Vc- This flux then fixes the flux Fj?' and the visibilities Vp of 
T Tau Sa. The final, best estimates of the fluxes for T Tau Sa and 
T Tau Sb are given in Figure 13 and the visibilities for T Tau Sa 
are displayed in Figure 14. 

Since the determination of the visibilities of T Tau Sa is cru- 
cial, we checked how much the result is dependent on the model, 
i.e., the visibilities, of T Tau Sb. Even when T Tau Sb would be 
unresolved by the VLTI, i.e., Vc = 1 -0, the derived visibilities Vp 
for T Tau Sa are almost indistinguishable from those that have 
been derived on the basis of the model for T Tau Sb. The same 
is true, when assuming that Vc is 0.6 and thus close to the mea- 
sured visibilities of the southern binary. A comparison is made 
in Figure 14. 
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Table 5. The parameters of the radiative transfer model of T Tau Sa. 
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6.3. A Radiative Transfer Model for T Tau Sa 

Following the approach of Duchene et al. (2005), we assumed 
that T Tau Sa has an effective temperature of 9000 K and a lumi- 
nosity of 40 L Q . The photosphere is represented in our model by 
a black body. At 2.3 M Q , T Tau Sa is the most massive star in the 
triple system (Kohler 2008). 

Similar to our approach for T Tau Sb, we modelled T Tau Sa 
as a source that is surrounded by a small truncated disc with a 
size of 5 AU (Table 5). This circumstellar structure contains a 
mass of 0.003 M Q of gas and dust and thus very efficiently ab- 
sorbs the photons from the central star due to the high inclination 
of the disc of around 70°. In addition, an envelope with a steep 
density gradient seems to be present. On the other hand, the con- 
tribution of accretion is negligible. Is has to be mentioned that 
without the foreground extintintion of 30 mag in the visual the 
silicate feature would again appear in emission like it has been 
found for T Tau Sb. 

In Section 6.1 for the two southern sources a preliminary 
optical depth of Ts a =1.8 and Tsb=0.9 was found. When fit- 
ting the spectra plotted in Figure 13 using the optical depths 
for astronomical silicates towards the Galactic centre (Kemper 
et al. 2004), one finds Ts a =l-7 and Tsb=0.8. Conversion of the 
silicate absorption depth into visual extinction, e.g., Rieke & 
Lebofsky (1985), leads to A v ~ 30 mag towards T Tau Sa and 
A v ~ 15 mag towards T Tau Sb. These results are consistent 
with the values used in our models and are attributed to the fore- 
ground screen. 

A possible explanation for the additional 15 mag of visual 
extinction towards T Tau Sa might be a circumbinary disc cur- 
rently only extincting this source. Such a scenario has been 
tested by reducing the foreground extinction towards T Tau Sa 
to 15 mag, extending its circumstellar disc to 50 AU, and cutting 



\ [nm] 

Fig. 15. The spectral energy distribution for the model of T Tau Sa with 
a disc inclination of 72° when neglecting the accretion (solid), removing 
the "envelope" (dotted), changing the outer radius of the disc to 4 AU 
(dashed), to 6AU (dashed-dotted), and using/? = 1.25 (dashed-dotted- 
dotted). 



off therein the area between 5 AU and 30 AU. In this model the 
whole circumstellar structure contained 0.02 M of gas and dust. 
It sucessfully reproduced the SED and the visibilities. 

The SED derived from our model of a compact circumstellar 
disc is shown together with the photometric data in Figure 14. 
Also the upper and lower limits of the modelled visibilities of 
both T Tau Sa and T Tau Sb are plotted. The visibilities for 
T Tau Sb were used in Section 6.2 to determine the individual 
spectra of the two components of the southern binary and to sep- 
arate the measured visibilities of T Tau S. The derived visibilities 
for T Tau Sa are indicated in Figure 14 by the error bars and fall 
right within the range predicted by the radiative transfer model. 
The two-dimensional N-band visibilities are shown in Figure 18 
and will be analysed later to determine the orientation of the cir- 
cumstellar disc (Section 7.3). 

In Figure 15 we show the spectral energy distribution when 
changing different parameters of the model while keeping the 
inclination fixed. This demonstrates that the model is to a cer- 
tain level degenerate. For example adjusting the inclination of 
the models with a different outer radius leads to similar good 
agreement with the SED than our best-fit model. Neither the N- 
band, nor the K-band visibilities can help to overcome this. On 
the other hand, removing the envelope and especially changing 
its radial density distribution (not shown) leads to unacceptable 
SEDs that cannot reproduce the N-band fluxes well. Also chang- 
ing p to a value of 1.25 is not reasonable. 

Finally, we report on the predictions made by our model 
for T Tau Sa for the K-band visibilities at 2.2/im. We find for 
the 85 m baseline used by Akeson et al. (2002) visibilities be- 
tween 0.03 and 0.38, while the visibilities for the 105 m base- 
lines lie between 0.01 and 0.29. The exact values depend on the 
actual position angles of the baselines. In Figure 16 the two- 
dimensional K-band visibility is presented. 

7. The T Tau system 

7. 1 . Three-Dimensional Geometry 

Figure 17 portrays the T Tau triple system schematically, but to 
scale. The graph is centred on T Tau N. 

The orbit of the centre-of-mass of T Tau S around T Tau N 
is drawn in grey and the positions of the centre-of-mass of the 
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Fig. 16. The K-band visibility of T Tau Sa at 2.2 //m as predicted by our 
radiative transfer model. 

southern binary in intervals of 50yrs are indicated by crosses. 
At the time of the interferometric observations with MIDI, the 
orbit of the two southern components T Tau Sa and T Tau Sb 
around each other is overlaid (see Figure 9). The positions of 
the two southern components are also shown for the time of our 
measurements. Both orbits are taken from Kohler (2008). While 
the southern orbit has a period of 28 [+12, -4] yrs, a semi-major 
axis of 13[+6,-2]AU, an eccentrity of 0.47 [+0.16, -0.18], 
and an inclination of 34°8 [+12°0, -31°0], the wide orbit has 
a period of 17500 [+7500, -16000] yrs, a semi-major axis of 




-1 O l .... i .... I .... i .... I 
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Fig. 17. A sketch of the T Tau system combining results of our MIDI 
measurements and previous observations. The dashed ellipse labelled 
'C indicates the emission structure found by Herbst et al. (2007). The 
orbits are taken from Kohler (2008). The dotted line represents the pe- 
riastron, the dashed-dotted line the line of nodes. See text for further 
details. 



1200 [+1200, -600] AU, an eccentricity of 0.92 [+0.07, -0.08], 
and an inclination of 56° 1 [+3°0, -4°0]. The wide orbit is much 
less constrained due to the fact that only a small part of it is 
covered by the near-infrared observations. 

The circumstellar discs around the three components of the 
T Tau system appear in black. The sizes and the inclinations of 
the discs are those that we derived from our simultaneous fits 
to the spectral energy distributions and the visibilities. The size 
of the disc around T Tau N thus might be similar to the perias- 
tron distance of the wide orbit. This indicates that the south- 
ern binary has to pass its periastron in the background or in 
the foreground with respect to the plane of the almost face-on 
disc around T Tau N. Otherwise dynamics would lead to a much 
smaller truncated disc. Surprisingly, the orbital elements of the 
wide orbit show that the southern binary is close to the line of 
nodes. The three stars thus nowadays reside at similar distances 
from the sun and the distance of the southern binary from the 
face-on disc has to be indeed quite small. Here it has to be men- 
tioned again that neither the wide orbit, nor the size of the disc 
around T Tau N is well constrained yet. So we have to leave it 
open, whether refined orbital elements, a smaller outer disc ra- 
dius or even a larger disc with a gap will help to overcome this 
inconsistency. 

Our simulations show that the two discs of the southern bi- 
nary are not aligned. While the inclination of the truncated disc 
around T Tau Sb is consistent with the inclination of the orbit 
of this source around T Tau Sa (Kohler 2008), the disc around 
T Tau Sa is highly inclined. The inclination of the disc around 
T Tau Sa is again in good agreement with the inclination of the 
wide orbit. Finally, the face-on disc around T Tau N exhibits a 
third orientation in the T Tau triple system. 

7.2. Formation of the system 

When fragmentation of a rotating cloud is responsible for the 
formation of a binary, the circumstellar discs are preferentially 
co-aligned with respect to each other. Also the orientation of the 
plane of the orbit should be similar to that of the discs. Therefore, 
one would intuitively argue that the misaligned discs in the T Tau 
triple system favour an independent formation of T Tau N and 
T Tau S. Consequently, even the misaligned two southern com- 
ponents might have formed independently. However, only in the 
centre of very dense clusters the required subsequent captur- 
ing process is efficient enough (Hall et al. 1996; Bonnell 2001). 
Formation by capture is thus very unlikely for T Tau. 

In order to identify a viable formation scenario we have to 
determine the most fundamental parameters of the system. One 
finds that about half of the mass of the system is concentrated in 
the northern, the other half in the southern component. In addi- 
tion, these two mass concentrations orbit around their centre-of- 
mass in a highly eccentric orbit. 

A scenario that leads to such a configuration is that of a fil- 
amentary collapse (Zinnecker 1989, 1991; Bonnell & Bastien 
1992). A gravitationally unstable, elongated fragment tumbling 
"end over end" will subsequently break up into two or even 
more fragments. Initially, these fragments are heading towards 
each other due to their gravitational attraction. But for a non- 
zero impact parameter the angular momentum will prevent the 
fragments from merging. Already formed discs will evolve and 
either directly or tidally interact with each other. 

The details of the process and the final configuration of the 
system are very sensitive to the initial conditions (Bonnell & 
Bastien 1992). For example, a wide range of binary mass ra- 
tios can be found in simulations of filamentary fragmentation. 



18 



Th. Ratzka et al.: Spatially resolved mid-infrared observations of the triple system T Tauri 



But the mass ratio tends to equalise during the evolution of the 
system, because the accretion onto the secondary is very effi- 
cient (Bonnell & Bate 1994). The secondary is able to collect 
even the high specific angular momentum material, as it resides 
in a circumstellar structure around the primary, away from the 
centre-of-mass. Furthermore, the eccentric orbit of the binary is 
circularised with time, because the specific angular momentum 
of the material accreted by the companion varies with the dis- 
tance from the primary. 

Another scenario to explain the actual properties of the T Tau 
system might be a three-body interaction within a flattened par- 
ential envelope (Reipurth 2000). In a close triple approach a non- 
hierarchical triple is tranformed into a tight binary and a third 
component on a long-period or even hyperbolic orbit. Such in- 
teractions can also form giant bow-shocks and a circumbinary 
disc is finally surrounding the binary. Also the IRC phenomenon 
can be explained by this scenario, when assuming that only the 
ejected companion has left the large flattened envelope while the 
binary is still hidden therein. 

Whether the T Tau triple system was really formed by fil- 
amentary fragmentation or a three-body interaction is unclear. 
A detailed investigation is far beyond the scope of this paper. 
However, only a highly dynamical process can be considered 
as a reasonable formation scenario (Pringle 1989; Whitworth 
2001). 

7.3. Origin of the jets 

To cite Bohm & Solf (1994) on the T Tau system: "It is intrigu- 
ing to see two distinct bipolar outflows of different orientation 
in the same stellar system". The jet pointing west and ending in 
the Herbig-Haro object HH 155 as its working surface (Biihrke 
et al. 1986), was attributed by them to T Tau N. This association 
appears convincing: The inclination of the jet axis against the 
viewing direction was determined from a comparison of radial 
and tangential velocities to a moderate value of 23° (Eisloffel & 
Mundt 1998) and could be seen in the forbidden lines of [OI] 
and [SII] well collimated to within 0.3" of the star. The star 
was found to be also mildly inclined by 19° (Herbst et al. 1997). 
Millimetre interferometry of its circumstellar disc found an incli- 
nation of i = 41° + 3° and a PA of 19° ±5° (Akesonetal. 1998). 
Near-infrared interferometry (Akeson et al. 2002) confirmed the 
inclination with i = 29° + 12° but disagreed on the PA with 132° 
+ 16°. Our model above, fitting SED and mid-infrared visibil- 
ities, led also to a moderate inclination of ; < 30°, while the 
PA could not be determined. With this evidence it was natural 
that Bohm & Solf (1994) attributed the second outflow, oriented 
nearly north-south, to the then known second component in the 
system, T Tau S. 

A Fabry-Perot map of the T Tau system in the v = 1 - 
0S(1) line of molecular hydrogen led Herbst et al. (2007) 
to contradicting conclusions. They found about one arcsecond 
west of T Tau S a bow shock-like emission, which the authors 
named "C". Their map clearly shows a huge complex consist- 
ing of a bow shock and the walls of the tunnel produced by the 
jet. This complex is indicated in Figure 17 by the dashed el- 
lipse. With the high spatial resolution provided by the instrument 
NACO on the VLT the jet could be unambiguously traced back 
to T Tau S, with the dominating component T Tau Sa being the 
preferred candidate for its origin. The northwest-southeastern 
(actually NNW-SSE) outflow, which is the inner part of the 
north-south oriented outflow mentioned above, could be best 
traced back to T Tau N. 



To increase the difficulty in combining these two sets of mea- 
surements and interpretations, Solf & Bohm (1999) found a high 
inclination of i = 79° for the north-south oriented outflow on the 
basis of the high dispersion/average value for its velocity pattern, 
which makes it hard to imagine the almost face-on circumstel- 
lar disc of T Tau N as its source. These complicated geometrical 
relations led Herbst et al. (2007) to list four different reasonable 
combinations of outflows with their stellar sources. 

To contribute to this discussion, we also have to derive the 
position angle of the disc around T Tau Sa in addition to the in- 
clination of 72° resulting from the simultaneous fit to the SED 
and visibilities. To this end, the modelled, two-dimensional visi- 
bilities for T Tau Sa are drawn in Figure 18. Since the orientation 
of a baseline with respect to the model, i.e., the object in the sky, 
is uncertain, the spatial frequencies to which a baseline is in prin- 
ciple sensitive fall onto a circle, whose radius is determined by 
the projected length of the baseline. However, the relative orien- 
tations of the three used baselines are fixed. Or in other words, 
the three baselines represented in the figure by lines from the 
origin to the circles can be rotated around the origin, but only 
together. On the other hand, only at certain locations along the 
circles the visibilty values are consistent with the visibilities de- 
rived from the measurements. These locations are marked by the 
thick white arcs in Figure 18. If it is now possible to find a rota- 
tion angle that brings all three baselines into marked locations on 
the corresponding circles, the orientation of the u-v-plane with 
respect to the baselines and thus to the sky is known with an 
ambiguity of 180°. 

The modelled visibility for T Tau Sa at a wavelength of 
8.4yLan is drawn in Figure 19 also as function of the position 
angle of the disc's minor axis. The model is calculated in such a 
way that without a rotation the minor axis of the disc is oriented 
north-south, and thus the two-dimensional visibility is elongated 
along the north-south direction. Therefore, the visibility predic- 
tion for the baseline with the longest projected length, which 
had a position angle of about 90° on the sky, has a minimum 
without a rotation of the model disc, and a maximum when the 
model disc is rotated to 90°. Similarly, the maxima for the other 
baselines also are reached when the model disc is rotated by an 
angle that corresponds to the position angle of that baseline. The 
visibilities derived from the measurements are drawn as bars in 
the right panel, and the resulting allowed ranges of position an- 
gels are marked in the graph by shaded areas (Figure 19). The 
height of such an area reflects the uncertainty of the correspond- 
ing measured visibility. 

It appears that agreement between the modelled and the mea- 
sured visibilities is found when the model disc is rotated by 
about 100°, i.e., when the disc is elongated almost along the 
north-south direction. This finding is robust, because the mod- 
elled visibilities for the different baselines in Figure 19 show 
slopes of opposite sign due to the phase shifts of the maxima. 
This may qualify T Tau Sa as a launching site for the east- 
western jet (Figure 17). Especially, the position angle of the 
bow-shock CI at the western tip of complex C (Herbst et al. 
2007) is in excellent agreement with the found position angle of 
the circumstellar disc around T Tau Sa. Also the above discussed 
scenario of a circumbinary disc only extincting T Tau Sa could 
be realised, because T Tau Sb would reside west of the north- 
south oriented disc plane. Interestingly, if such a circumbinary 
disc is present, T Tau Sb would be regularly extincted by this 
structure due to its orbital motion. 

If the determined north-south orientation of the circumstellar 
disc around T Tau Sa is correct, T Tau N would be the origin of 
the northwest-southeastern outflow with the usual reasoning by 
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Fig. 18. Top: The modelled visibilities of T Tau Sa for 8.4yum, 10.6/im, and 12.6//m and a disc inclination of 72° (from left to right). The solid 
contours are drawn for steps of 0. 1 and the associated grey background is underlaid for easier comparison. The radii of the overplotted circles 
correspond to the spatial frequencies for which the baselines are sensitive at the various wavelengths (solid: 88 m projected baseline length, dashed: 
62 m, dotted: 44 m). Bottom: Zoom in on the panels in the top row and determination of the orientation of the u-v-plane relative to the sky (see 
text). 



exclusion. Unfortunately, the difficulties with picturing the T Tau 
system of flows have not been reduced. The found orientation for 
the disc of T Tau Sa indeed perfectly matches the east-western 
flow of Figure 17 if it is close to the plane of the sky. But then the 
difficulty arises to combine this flow with the only moderately 
inclined jet seen at larger distances out to HH 155. The same 
inclination discrepancy emerges for the north-south going out- 
flow. To postulate that orientations may have changed strongly 
in this triple system during the about one million years it may 
have taken the jet (v ra d = 39km/s) to arrive at the position of 
HH 155 has no serious foundation at present. On the other hand, 
the length of a Kozai cycle for the secular evolution of the incli- 
nation of the southern binary should be comparable (Ford et al. 
2000) and episodic ejections are much shorter (Reipurth 2000). 
Perhaps a kinematical determination of the outflow inclination 
from the H2 emission west of T Tau would be a first clarifying 
step. Further questions would concern better knowledge of the 
three-dimensional grouping of outflows around T Tau and a bet- 
ter understanding of the role of the "second" infrared companion 
T Tau Sb. There is no reason why T Tau Sb with its circumstellar 
disc and the highest accretion rate among the three stars in the 
T Tau system should not be the origin of one of the two jets. 

We presented new information for discussing this puzzle of 
the origin of the jets, a solution still has to be awaited. 

8. Summary 

MIDI at the VLTI allows an interferometric study of the three 
individual components of the T Tau system and their associated 
circumstellar material. 
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Fig. 19. The modelled visibilities of T Tau Sa drawn as a function of the 
position angle of the minor axis of the disc (solid: 44 m projected base- 
line length, dashed: 62 m, dotted: 88 m). For comparison, the expected 
visibilities are plotted in the right panel. The corresponding shaded ar- 
eas in the graph indicate the suitable position angles. 



The disc around T Tau N was well resolved in our obser- 
vations. We find a radiative transfer model, including accretion, 
that fits simultaneously the SED and the visibilities. The physi- 
cal parameters are in line with earlier measurements. This means 
that we can confirm the picture of an almost face-on, compara- 
tively massive circumstellar disc around T Tau N. An optically 
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thin envelope is needed as an additional component to produce 
the 3-8 yum radiation. The model also fulfils the constraints set 
by the K-band visibilities measured with the PTI (Akeson et al. 
2002). By fitting the silicate emission feature in the total spec- 
trum and the correlated spectra, we find that dust processing oc- 
curs in the inner parts of the circumstellar disc around T Tau N. 

The binary T Tau S is also well resolved interferometrically, 
but not photometrically. The visibilities with the imprinted si- 
nusoidal modulation give the relative position of T Tau Sa and 
T Tau Sb and allow us to derive the flux ratio of the (corre- 
lated) fluxes. The silicate band is seen in absorption towards 
both sources and confirms a high foreground extinction of about 
A v = 15 mag. The phases indicate that T Tau Sa is the brighter 
source in the N-band. 

For the full use of the data we adopt for T Tau Sb the model 
of Duchene et al. (2005), which describes this source as a normal 
T Tauri star of spectral type Ml. The star resides behind the ab- 
sorbing screen of Ay = 15 mag and its circumstellar disc is seen 
not far from face-on. The low-mass disc is truncated because of 
the proximity of T Tau Sa. 

With this model at hand, the interferometric data reveal an 
elongated structure for T Tau Sa, interpreted as a disc seen nearly 
edge-on. The radiative transfer model with a central star as pro- 
posed by Duchene et al. (2005) simultaneously fits SED and vis- 
ibilities. The model again suggests the presence of an optically 
thin envelope. The orientation of the disc is almost north-south. 
T Tau Sa might thus be the driving source of the east-west jet. 
The additional (intrinsic) extinction of about A y = 15 mag might 
be attributed to a circumbinary disc. 

The orbits and discs in the T Tau triple system show complex 
orientations. While the disc around T Tau Sb might be coplanar 
with its orbit around T Tau Sa, the disc around T Tau Sa is not 
aligned. It is almost perpendicular to the disc around T Tau N. 
However, the full picture requires knowledge of the orbital ele- 
ments of the orbit of the southern binary around T Tau N with 
much better precision than is available now. 



Appendix A: Visibility of a binary with extended 
components 

Single-baseline interferometry gives one-dimensional informa- 
tion on the object along the direction of the baseline projected 
onto the sky. Let P A and C A be the spectrum of the primary and 
the companion of the binary with projected brightness distribu- 
tions along this direction of P,i(x) and C A (x). The complex vis- 
ibilities, the Fourier transforms with respect to x, are written as 
P A (u) and C A (u). The correlated fluxes of the components then 
are F c p 0I (u) = \P A (u)\ and F c c or (u) = \C A (u)\. The normalised 
visibilities, also functions of A, are V p (u) = \P A (u)\ jP A and 
V c (u) = \C A (u)\ IC\. Please note that in the case of a measure- 
ment with MIDI, u and A are coupled by the relation u = B/A, 
where B is the length of the projected baseline. For extraction of 
the phase from MIDI measurements, see Appendix B. 

In the special case where the baseline projection is perpen- 
dicular to the separation vector of the binary, we have 

1. the flux distribution A (x) = P A (x) + C A (x) , 

2. the correlated flux F c ° r = \P A (u) + C A (u)\ , and 

3. the normalised visibility V„ = ^ffig 001 ■ 

If P A (x) and C A (x) are symmetric (Fourier transform is real) 
or the same apart from a constant factor, or if P A (u) and C A (u) 



have the same phase, then we have the naively expected result 
that the correlated flux of the binary is the sum of the correlated 
fluxes of the components 



\r A { U )\ + \c A {u) 



(A.l) 



and the normalised visibility of the binary is the flux weighted 
average of the visibilities of the components 



V («) 



p A y p (u) + c A v c (u) 
p a + c a 



(A.2) 



It is common to the three conditions given above that they pre- 
sume the same phase for V A {u) and C A (u). If this does not hold, 
relations (A.l) and (A.2) are not fulfilled in the strict mathemat- 
ical sense. But they still may serve as good approximation, to be 
judged on a case-by-case basis. 9 

For all other projected baseline directions, the primary and 
companion are separated by a projected separation s. The corre- 
lated flux is then 



f 



[P A (x) + C A (x - s)] e- lniux dx 



or after substituting x-s 

F c m (u) = \<P A (u) + e- 2 " ius C A (u)\ 



(A.3) 



(A.4) 



= y/(P A (u) + e- 2nius C A (u))(P A (u) + e- 2nius C A {u))* , 

where the asterisk denotes the complex conjugate of the expres- 
sion in parentheses. 

Writing <P A (u) = \V A (uy\e^' (u) and C A (u) = \C A (u)\ e* M , and 
with A<f> A (u) = </> p (u) - <f> c (u), the argument of the square root has 
the value 

\P A (u)\ 2 + \PAu)\ \C A {u)\ cos(AMu) + 2ms) + \C A {u)\ 2 . (A.5) 

For s — this is the detailed formulation for the correlated flux 
observed perpendicular to the separation vector of the binary. For 
A(p(u) = and constant values of V A (u) and C A (u), this reduces to 
the well-known formula for a binary composed of point sources. 
As in the case of point sources, in this more general case the 
cosine term is responsible for the cosine-like modulation of the 
visibility, the correlated flux, respectively. 

It is a matter of convenience whether one analyses the corre- 
lated flux with the equations given here or the normalised visibil- 
ity. Then the given equations have to be divided by the combined 
flux of the binary. In the case of T Tau S, where the visibility 
varies smoothly with spatial frequency u and hence A, while flux 
and correlated flux show a strong depression in the silicate fea- 
ture, it is advantageous to discuss and interpret the visibilities. 

In the following paragraphs we will now consider two spe- 
cial cases for the product u ■ s. 



An integer result for the product u ■ s The cosine factor equals 1 
for A<f> A (u) = 0. This corresponds closely to the maxima of the 
modulated visibility curve. Here, the correlated flux and hence 
the visibility are exactly the same as for observations perpen- 
dicular to the separation vector of the binary. For A<p A (u) = 0, 
or if A<p A (u) + 2nus is an integer multiple of 27r, the simple re- 
lations (A.l) and (A.2) hold as given above. And, at least for 
A<p{u) = 0, the gradient in these maxima is the same as for the 



9 In fact, we will see that for a given wavelength, the phase equality 
necessary for strict validity of relations (A. 1) and (A.2) will be fulfilled 
at suitable values of the spatial frequencies. 
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unmodulated curve that would be obtained perpendicular to the 
separation vector. Therefore, to the extent that linear interpola- 
tion of the visibility appears justified, the combined visibility at 
the given baseline can be reconstructed as a function of wave- 
length from the two maxima appearing in the modulated visibil- 
ity curves observed for T Tau S. 

Case of a half-integer result for the product u-s The cosine 
factor equals -1 if A(f> A (u) + 2nus is an odd multiple of n. This 
corresponds closely to the minima of the modulated visibility 
curve, with F z °\u) = \P A (u)\ - \C A (u)\ for A<f> A (u) = 0. 

Comparing this minimum of the correlated flux or visibil- 
ity to the maximum value interpolated to the same spatial fre- 
quency for the unmodulated visibility curve, we have essentially 
the same relation as for simple binary stars: 



F c 0I (max) - F™ r (min) 
F™ r (max) + F c or (mm) 



V (max) - V (min) 
V (max) + V (min) 
= 2 ■ CM 
' 2-<P A (u) ' 



(A.6) 



This means that to a good approximation, we can read from 
the observed modulated visibility curves the ratio of the corre- 
lated fluxes of the components at this spatial frequency. By fit- 
ting the shape of the observed visibility modulation for the en- 
tire available wavelength range, this ratio can be inferred as a 
function of u. In binaries composed of point sources, the terms 
in which we are used to think, the equivalent relation gives the 
brightness ratio of the components, because flux and correlated 
flux are the same for point sources. With the assumption that the 
primary and companion show the same visibility, V p (u) = V c (u), 
Equation (A.6) would also give the brightness ratio. 



Appendix B: Appendix B. Phase measurements 
with MIDI 

The basic equations Consider an object P A (a), where a is an 
angle on the sky. We write its one-dimensional Fourier trans- 
form, i.e., its complex visibility, as 

PAa)e- 2mua da 

oo 

= \<P x (u)\ ■ e* (u) = A A (u) • e* (u) 



(B.l) 



with visibility amplitude A A (u) and phase <p A {u). Reversing the 
direction in which a is increasing is equivalent to going from 
P A (a) to P A (-a) and would change the sign of <p A {u). Observing 
this object interferometrically with a baseline of projected length 
fiproj and a corresponding spatial frequency u = B pm j/A on the 
sky leads for multiaxial beam combination to a visualisation of 
the visibility component P A (u) as an instantaneous fringe pattern 
in the image plane with the modulating part 



h(x) oc 2\P A \cos(2jrf m x + 4> A ), 



(B.2) 



where l// m is the spatial fringe spacing and the sign of the phase 
is positive for proper choice of the sign of x. Alternatively, in 
coaxial beam combination, this visibility component shows an 
OPD-scanned timewise fringe pattern of 



I A (OPD) oc 2|^|sin(27rOPD//l + i ), 



(B.3) 



VLTI. The factor of two here results from taking the difference 
of the two beamsplitter output signals of opposite sign. The pres- 
ence of the sine function in relation (B.3), instead of the cosine 
function, results from the well-known n/2 phase shift in the sym- 
metrically used beam splitter, which allows the coaxial recombi- 
nation. As we will see, the n/2 phase shift is of no consequence 
for the phase determination of <f> A . The positive directions of x 
in Equation (B.2) and OPD in Equation (B.3) need to corre- 
spond to the positive scan direction a over the source on the 
sky. Otherwise, the phase <f> A would change sign. In other words, 
defining the signs of x, OPD, respectively, and keeping Equation 
(B.2), (B.3), respectively, would also define a positive direction 
a on the sky. 

For later use we note that by the formulation of 
Equation (B.2) the maximum of the fringe pattern at the origin 
will move to negative values of x when cf> A is going positive from 
the value (which a point source at the origin would have). In 
the coaxial case the zero crossing of the fringe pattern at the ori- 
gin will, according to Equation (B.3), move to negative values of 
OPD for an increasing <p A . 

Derivation from more general principles Equations (B.2) and 
(B.3), though plausible, deserve some justification in detail, in 
particular, to verify the sign of the angle a, of the OPD, and of 
the phases. We start from the presentation by Boden (2000). A 
plane wave propagating from the object in direction k towards 
the observer is written as e lVjL ~ ,M with |k| = k = 2n/A. Using 
telescopes A and B we define the baseline vector B as A — > B. 
Pointing the telescopes to the object, we are viewing in direc- 
tion s = -k. The total optical path difference counted in the 
sense B minus A then is OPD = -sB + ds - d A , where d are 
the pathlengths from the telescopes to the locus of beam combi- 
nation and -sB constitutes the external part of the OPD, which 
decreases for viewing directions towards telescope B. From the 
superposition of the waves travelling through telescope A, i.e., 
e i(sB+d A )-icot^ and te i esc0 p e b, i.e., e KdB ~ ia "\ the fringe pattern is 



F^(OPD) oc cos [k (sB + d A - d B )] ■ 



(B.4) 



where obviously AOPD = A is the fringe spacing during an OPD 
scan. The latter method is used in the MIDI instrument at the 



To relate positions on the sky to the instrumentally available in- 
ternal OPD d A - dg, we consider the viewing direction s as the 
origin on the sky. We take as the zero point of OPD that value 
of dA - ds where a point source at this origin would produce 
the white light fringe, i.e., where the argument of the cosine 
in Equation (B.4) is zero. We also note that with the accepted 
sign convention for the OPD (OPD = path length "MIDI beam 
B" minus path length "MIDI beam A"), an interferometric scan 
with MIDI moves the internal OPD ds - dA stepwise towards 
more positive values, since the piezo motion is shortening d^. In 
the current MIDI mechanical-geometrical setup, during such an 
OPD scan the corresponding position of the white light fringe on 
the sky moves towards the telescope that sends its light beam into 
the right entrance window of the MIDI dewar, when seen from 
the viewpoint of the incoming light ("MIDI beam B"). The FITS 
headers of the stored data files carry the information that allows 
transposition of this somewhat convolved definition into a well- 
defined direction on the sky, such that the phases will give well- 
defined information on the object. It is natural to count a pos- 
itive on the sky in the direction where the corresponding white 
light fringe moves during an internal OPD scan. According to 
Equation (B.4), this means that, because d A decreases during the 
scan, a as sB goes positive towards telescope B. 

For an extended source, we have to add the independent 
fringe patterns (Equation B.4) contributed from all its parts, 



22 



Th. Ratzka et al.: Spatially resolved mid-infrared observations of the triple system T Tauri 



which means integrating over a, i.e., in the plane of the sky and 
perpendicular to s. The integrand then reads 



P A (a) ■ cos [k (s + a) B - k (d B -d A + AOPD int )] . 



(B.5) 



Here, ds - flU is the value of the internal OPD corresponding to 
the white light fringe at a = 0, and AOPDi nt is an instrumental 
offset from the white light position (describing, e.g., the piezo 
scan). Since sB + dU - d B = 0, the integrand simplifies to 

P A (a) ■ cos [k(aB - AOPD int )] = 

P A (a) cos(27rafl proj M) cos(27rAOPD int //i) + 

P A {a) sin(27raB pr0 j /A) sin(27rAOPD int /i) . (B.6) 

Keeping in mind the above definition of visibility P A (u), and 
respecting <f> A = 2naB pro j/A, the integral over a refers to the in- 
tegrand 

\T A (u)\ cos(4> A ) cos(27rAOPD int /T) - 

\<P A (u)\ sin(^) sin(2^AOPD int /i) , (B.7) 



which directly leads to the form 

^(AOPDi nt ) = \<P A (u)\cos(2nAOPD m / A + (/>,), 



(B.8) 



from which Equation (B.3) follows simply by applying the n/2 
phase shift. This verifies the formulation of Equation (B.3), in- 
cluding the sign of the phase. 

A Gedankenexperiment We now consider a binary with the 
centre of light in the origin to confirm by a detailed look at the 
expected temporal fringe pattern that the phases determined by 
MIDI will have the correct sign. 

On the sky we are on safe grounds: The phase would show 
the form of a staircase with width of the individual steps in spa- 
tial frequency u of Am = 1 / separation and the height of the steps 
depending on the brightness ratio. The phase will decrease with 
u if the companion is located at more negative a with respect to 
the primary. On the other hand, the binary's phase with removed 
linear trend - as it is determined in the instrument MIDI - cor- 
responds to a different position, namely with the primary at the 
origin. The binary phase now oscillates around zero. For a com- 
panion at more negative position a than the primary, the phase 
will increase with u (decrease with A) at those spatial frequency 
values where the maxima of the visibility amplitude are found, a 
simple and useful criterion. 

In the instrument it looks like this: If we have a companion at 
slightly negative a with respect to the primary, its contribution 
to the temporal fringe would have a maximum slightly earlier 
in the OPD scan, at AOPD comp = -a ■ B. The maximum of the 
combined fringe pattern (primary plus companion) then will also 
be shifted to more negative values, by 



AOPD binaiy = -a ■ B ■ 



comp 



^prim + ^comp 



(B.9) 



Measured in units of A, this shift, having a negative value, will 
get closer to zero, i.e., will be rising with A. The phase <f> A , corre- 
sponding to this shift of maximum, therefore will be decreasing 
with A. This relation has the correct sign, and so the phases deter- 
mined from MIDI should be correct with the sign as defined by 
the Fourier transform and as used in Equations (B.3) and (B.8). 



The treatment of phases in MIDI In MIDI, OPD in 
Equation (B.3) is the actual optical path difference for a given in- 
dividual exposure. One such exposure with the MIDI instrument 
gives /^(OPD) as function of A for the range 8-13/im. Many 
of these exposures at a repeated saw-tooth pattern of OPD steps 
scanning the fringe constitute an interferometric measurement. 
In each scan, the position of fringe maximum depends on the 
phase <f> A . 

For the determination of phases from MIDI measurements 
we have to recall that the total optical path difference consists of 
two parts of different behaviour: OPD tota i = OPD; nt + OPD atm , 
where the "instrumental" OPD; nt is quickly stepped through a 
few times the wavelength of about 10 //m during an interferomet- 
ric measurement and is known including its sign, while OPD atm 
is mostly determined by the atmosphere and is more slowly vary- 
ing with a few /j.m/s, but in an unpredictable way. 

The EWS part of the MIDI data reduction package 
MIA+EWS (Section 3) determines the phase by mathematically 
shifting each exposure as given by Equation (B.3) to OPD zero 
and then averaging to obtain the complex visibility - modu- 
lus and phase (Jaffe 2004). Remembering that a shift in image 
space (and here also in OPD) corresponds to a phase factor in 
Fourier space of the form e 2mu0PD ^ we multiply (B.3) by the fac- 
to "shift" each exposure F A (OPD) back to posi- 



tor e 



-2m'MOPD in 



tion OPD; nt = 0. Equation (B.3) now reads - expressing the sine 
function by exponentials and apart from factors that do not mat- 
ter here - as follows: 



i^(OPD) = -j-e 



-1 „2m«OPD alm +(0j_2£ e -4ff(«OPD in ,-2^!«OPD a , m -/f ( jq) 



To determine the unkown atmospheric optical path differ- 
ence OPD atm , we Fourier transform Equation (B. 10) with respect 
to u and obtain peaks at +OPD atm and (-OPD atm - 20PDi nt ). Of 
these, the second term, rapidly changing with OPDi nt during a 
scan, quickly smears out by averaging over different exposures, 
so the only remaining peak gives the wanted value of OPD atm 
with the correct sign. We can now perform the "shifting back" 
correction also for the atmospheric part of the OPD by multiply- 
ing Equation (B. 10) with e -^'"OPD atm This leads tQ 



/^'(OPD-corrected) = A A e +i( ^ ±n/2 \ 



(B.ll) 



where we have omitted the strongly jittering part with the factor 
g-toOPDint-toOPD^-* The constant results from the factor 

1/z. Its sign depends on which of the two complementary inter- 
ferometric signals after the beam splitter is subtracted from the 
other one. Averaging Equation (B.ll) over the typically several 
thousand exposures of an interferometric measurement thus will 
reliably determine - apart from a constant known in principle 
- the desired phase <f> A (u = B pro j/A) over the wavelength range 
8-13 iim U) . 

We remove any linear trend with u from this phase, since it 
would simply correspond to a position offset of the source from 
the origin. In addition, we remove any remaining constant term 
+0. Again, apart from our constant, this term could be and ap- 
pears to be mostly due to differential longitudinal dispersion over 
the unequally long air paths of the interfering light beams and 
it thus mostly carries information not related to the object. Of 
course, subtraction of the constant term also takes care of the 
added + n/2. This ends the determination of the phase. We note 
that the data reduction procedure described here preserves the 
sign of the phase as defined by the Fourier transform and given 
in Equation (B.3). 



At the same time, it will also give the visibility amplitude A A 
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The phases determined with MIDI on the VLTI thus are not 
fully complete, i.e., the constant and linear terms are set to zero. 
But they allow to resolve the ambiguity of where the companion 
in a binary is located with respect to its primary. This property of 
the MIDI data is needed and used in the discussion of the close 
binary T Tau S (Section 6.2). 

A check by a known binary We wanted to check the signs of the 
phase derived above from "first principles" with an object where 
we assume to know the position of the companion. The close 
binary Z CMa, with separation of 0.1", consists of an infrared 
companion and a FU Ori component (Koresko et al. 1991). At 
4 /an the infrared companion is brighter by a factor of 6. This ra- 
tio rises with wavelength. The ratio of correlated fluxes at lOyum, 
determined as described in Appendix A, is in the range 5-10. We 
assumed that with these clear inequalities the component much 
brighter at 4 /mi will also be the much brighter component at 
lOyum and that it also will show the larger correlated flux. Z CMa 
then can be used for calibrating the sign of the phase relation, 
with the result that the relations derived in this appendix appear 
correct. A final confirmation came from observations of r\ Vir 
(Appendix C). 

Appendix C: 'The r\ Vir Experiment' 

We got permission from ESO to prove experimentally the con- 
clusions derived in Appendix B. As target we chose the well- 
studied triple system 77 Vir (Hartkopf et al. 1992; Hummel et al. 
2003). T] Vir has a composite spectral type A2IV and is formed 
by a close spectroscopic binary and a wide companion. The or- 
bit of the spectroscopic binary has a semi-major axis of 7 mas 
and can thus only be marginally resolved by means of optical 
long baseline interferometry in the mid-infrared. On the other 
hand, the wide companion is in an orbit with a semi-major axis 
of 1 34 mas and can be well detected with a suitable baseline con- 
figuration. 

The observations were performed on May 17, 2006 with the 
baseline UT3-UT4. The projected length was 53 m at a position 
angle of 113°. According to the orbital parameters published by 
Hummel et al. (2003), the wide companion was at that time at a 
position angle of 28 1 and separated by 89 mas from the primary. 
The projected separation along the projected baseline was thus 
87 mas. For the calibration of the data, observations of e Crt have 
been obtained right after those of r\ Vir. The projected baseline 
length was 45 m at a position angle of 130°. 

In Figure C.l, the normalised visibility (solid black) and the 
phase (dashed black) are shown. The modulation of the visibil- 
ity is obvious, although the extremes are shifted to higher spatial 
frequencies when compared to the predictions by the orbital fit 
(Hummel et al. 2003). The latter are indicated by arrows. The 
MIDI measurements are best represented by a binary with a sep- 
aration of about 80 mas and a flux ratio around 0.25. The vis- 
ibility of such a binary model is shown as solid grey line in 
Figure (C.l) after multiplying it by a value of 0.8. The phase, 
plotted as dashed black line, shows a rising plateau (dashed grey) 
that is centred around the visibility maximum. The smooth tran- 
sition to the next plateau occurs as expected around the visibility 
minimum. The change in phase can be well estimated, although 
the second plateau can be seen only marginally. We measure a 
value of Aip « 75°. On the other hand, the value derived from 
the flux ratio of 0.25 is Acf> = 2n - 2n/(l + 0.25) = 2/5tt = 72°, 
almost identical. We thus conclude that we see in Figure C. 1 the 
signal of the wide binary. 
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Fig. C.l. The normalised visibility (solid black) and the phase (dashed 
black) of 7] Vir. The arrows indicate the predicted positions for the max- 
imum and minimum of the visibility, while the grey solid line represents 
a binary model. The grey dashed line indicates the plateau of the phase. 

With this knowledge, the conclusion drawn in Appendix B 
can now be tested. From Figure C.l, one sees that the phase 
increases with spatial frequency (decreases with wavelength) 
at the visibility maximum. According to Appendix B, this is a 
clear sign that the companion was scanned before the primary. 
Looking at the headers of the files one finds that beam A was 
associated with UT3, while beam B is related to UT4. Since the 
scan is done towards beam B, its direction was towards posi- 
tion angle 113°. This leaves us with the result that the compan- 
ion was indeed scanned first and that the conclusion drawn in 
Appendix B is correct. 
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